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ABSTRACT 

Optical  beam  pointing  is  a  critical  topic  in  the  study  of  Directed  Energy  Weapons 
Systems.  One  of  the  main  operational  concerns  with  optical  beam  pointing  is  the  effect  of  small 
vibratory  motion,  defined  as  jitter.  Understanding  jitter  is  important  for  minimizing  the  effect  it 
has  upon  a  directed  energy  beam’s  intensity  at  its  target.  In  this  project,  Lagrange’s  equations  of 
motion  are  derived  for  a  visco-elastic,  point- supported  plate  containing  discrete  masses,  which  is 
representative  of  an  optical  platform.  The  mechanical  response  of  the  optical  platform  is  sought 
for  free  vibration,  impact  loading  and  the  response  imparted  by  two  inertial  actuators.  All  three 
of  these  responses  are  then  analyzed  with  and  without  damping.  A  solution  of  these  equations 
contains  information  pertaining  to  the  plate’s  amplitude  and  frequency,  analyzed  in  both  time 
and  frequency  domains,  which  together  computes  the  effect  of  jitter  on  the  intensity  of  a  directed 
energy  beam  at  the  target. 


Keywords:  jitter,  Lagrange’s  equations,  viso-elastic,  actuators 
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NOMENCLATURE 

=  amplitude 
=  damping  coefficient 
=  critical  damping 

=  Elastic  modulus  (Young’s  modulus) 

=  cycles  per  second  (Hertz) 

=  mass  moment  of  inertia 
=  stiffness  coefficient 

=  mass 

=  Force  (Newton) 

=  matrix  of  eigenvectors 
=  Force  per  unit  area  (Pascal) 

=  Rayleigh’s  dissipation  function 
=  displacement  in  modal  space  with  respect  to  time 
=  transformation  matrix  from  modal  space  to  original  space 
=  time 

=  Period  of  oscillation 
=  Kinetic  Energy 
=  damping  ratio 
=  Spring  Energy 
=  damped  natural  frequency 
=  natural  frequency 
=  phase 
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CHAPTER  1 

INTRODUCTION  AND  BACKGROUND 

1.1  Introduction 

Directed  energy  systems  have  eluded  researchers  and  engineers  for  decades  due  to  the 
vast  difficulties  ranging  from  the  basic  physics  of  the  system  to  modeling  and  system  design. 
Advances  in  laser  and  control  systems  have  proved  fruitful  towards  solving  these  problems. 
Attention  is  now  turning  to  other  issues  that  affect  the  directed  energy  system.  One  such 
challenge  is  to  gain  a  better  understanding  of  induced  structural  vibrations  on  the  mechanical 
response  of  the  structure  that  contains  the  directed  energy  system.  This  vibration,  termed  jitter, 
has  a  measurable  effect  on  the  accuracy  and  intensity  of  the  laser  at  its  target.  A  model  capable 
of  simulating  the  motion  of  an  optical  platform  in  various  scenarios  would  greatly  assist  in  the 
understanding  of  jitter  effects  and  the  mechanics  behind  the  structural  support  system  of  the 
laser. 

1.2  Background 

The  optical  laser  system  can  be  modeled  as  a  multiple  degree  of  freedom  rigid  body  with 
attached  discrete  masses  supported  on  visco-elastic  springs.  One  approach  is  to  model  it  as  a 
flexible  body  and,  due  to  the  support  conditions,  include  rigid  body  modes.  Work  on  the 
vibration  of  deformable  plates  is  mature  and  only  an  overview  is  presented  here. 

Cox  and  Boxer  [1]  provided  one  of  the  early  contributions  to  this  topic  and  discussed  the 
vibration  of  a  point  supported  plate  at  its  corners  and  free  edge  boundaries  elsewhere.  Both 
rotary  inertia  and  shear  effects  were  neglected  for  the  plate.  A  finite  difference  approach  was 
used  to  solve  the  governing  differential  equation,  and  provided  mode  shapes  and  frequencies  for 
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both  square  and  rectangular  plates.  Amba-Rao  [2]  studied  the  vibration  of  a  simply-supported 
rectangular  plate  carrying  a  concentrated  mass  and  presents  a  closed-form  solution  for  the 
frequencies  and  mode  shapes.  To  develop  the  formulae  for  the  frequencies,  the  transverse 
displacement  was  represented  as  a  double  infinite  sine  series.  Kerstens  [3]  incorporated  a  modal 
constraint  method  to  solve  the  problem  of  vibration  of  a  rectangular  plate  supported  arbitrarily. 
The  modal  constraint  approach  used  the  principle  of  minimum  potential  energy,  which  included 
virtual  work  done  by  homogeneous  boundary  conditions.  The  point  supports  were  represented 
through  the  virtual  work  by  a  homogenous  boundary  condition  formed  by  finite  constraints  as 
point  supports.  A  Lagrange  multiplier  was  used  to  represent  the  unknown  forces.  The  use  of  the 
Rayleigh-Ritz  method  continues  to  be  a  common  approach  for  solving  the  eigenvalue  problem 
associated  with  plate  vibration.  Of  particular  interest  was  the  selection  of  basis  function  used  in 
the  analysis.  Lee  and  Lee  [4]  incorporated  a  new  class  of  admissible  functions  to  study  the 
vibration  of  elastically  point  supported  plates.  The  basis  function  corresponded  to  a  similarly 
supported  beam  under  point  loads.  Singularity  functions  were  used  to  completely  describe  the 
equation  defining  the  elastic  curve.  Normalized  frequencies  were  presented  for  a  plate  simply 
supported  on  all  sides  and  as  a  function  of  support  location  as  well  as  elastic  support  stiffness. 
Kim  and  Dickerson  [51  incorporate  orthonormal  polynomials  as  the  basis  functions  in  the 
Rayleigh-Ritz  method  to  study  the  vibration  of  point  supported  plates  with  a  combination  of 
simply-supported  and  clamped  edge  conditions.  The  orthronormal  set  was  constructed  using  the 
Gramm-Schmidt  process.  Kocaturk  et  al.  [6]  considered  the  dynamic  response  of  a  thin 
orthotropic  plate  supported  on  viscoelastic  supports,  which  all  had  stiffness  and  damping 
properties.  The  plate  was  subjected  to  harmonic  input  at  its  center  and  supported  added  masses. 
The  arrangement  of  the  masses  was  placed  symmetrically  along  the  diagonals  of  the  plate.  A 
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Ritz  method  was  used  to  compute  the  eigenparameters  presented  in  the  Lagrange’s  equation  of 
motion.  Results  were  presented  for  several  material  constituents  and  considered  symmetric 
modes  of  vibration  consistent  with  geometry.  The  authors  also  presented  results  for  the  force 
transmissibility  as  a  function  of  plate  material  properties  and  added  mass. 

Of  particular  interest  was  experimental  work  that  had  been  done  in  this  area  of 
investigation.  Although  not  as  voluminous  as  numerical  and  analytical  work,  several 
contributions  were  noteworthy  in  literature.  Nieves  [7]  et  al  used  laser  interferometry  to  identify 
the  flexural  vibrations  of  thick  isotropic  plate.  The  technique  was  said  to  consist  solely  of  out- 
of-plane  displacement  flexural  modes.  Analytical  results  from  the  Ritz  method  were  used  as 
comparisons  for  accuracy.  A  Mindlin  plate  theory  was  the  basis  for  the  analytical  study 
presented.  Lee  and  Kam  [8]  considered  the  effect  of  elastic-edge  supports  on  the  vibration  of 
composite  plates  as  a  means  of  determining  the  mechanical  properties  of  a  plate.  Impulse  data 
provided  the  necessary  frequencies  to  extract  the  mechanical  properties. 

Rigid  body  mechanics  offered  a  much  more  suitable  alternative  in  light  of  design, 
fabrication  and  construction  of  mounted  systems  aboard  ships  and  aircraft.  In  this  approach, 
kinematic  equations  predicted  motion  leading  to  an  understanding  of  the  vibration  response. 

1.3  Approach 

The  analysis  of  a  vibrating  system  usually  involved  mathematical  modeling,  derivation  of 
the  governing  equations,  solution  of  the  equations,  and  interpretation  of  the  results.  A  brief 
summary  of  each  step  follows.  More  detailed  descriptions  are  included  in  chapters  (2-6). 


Step  1 :  Mathematical  Modeling 
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The  purpose  of  mathematical  modeling  is  to  represent  all  the  important  features  of 
the  system  for  the  purpose  of  deriving  analytical  equations  governing  the  system’s 
behavior.  The  mathematical  model  may  be  linear  or  nonlinear. 

Step  2:  Derivation  of  Governing  Equations 

Several  approaches  are  commonly  used  to  derive  the  governing  equations  of 
motion.  Among  them  are  Newtonian  Mechanics,  Hamilton’s  Principle,  the  Principle  of 
Virtual  Work,  and  Lagrange’s  Equation.  When  only  single  degree  systems  are 
considered,  Newton’s  law  is  useful,  but  in  more  practical  problems  with  many  degrees  of 
freedom,  energy  considerations  can  be  combined  with  the  concepts  of  virtual  work  to 
produce  Lagrange’s  equations.  Lagrange’s  equations  provide  an  energy-based  alternative 
to  summing  forces  to  derive  equations  of  motion  and  are  therefore  used  in  this  project. 

Step  3:  Solution  of  the  Governing  Equations 

The  equations  of  motion  must  be  solved  to  find  the  response  of  the  vibrating 
system.  Depending  on  the  nature  of  the  problem,  one  of  the  following  techniques  may  be 
used  to  find  the  solution:  standard  methods  of  solving  differential  equations,  Laplace 
transform  methods,  matrix  methods,  and  numerical  methods.  The  resulting  problem  for 
the  vibration  of  the  rigid  problem  is  an  eigenvalue  problem.  Modal  analysis  is  one 
approach  to  solve  the  associated  eigenvalue  problem  using  vector  transformations. 

Step  4:  Interpretation  of  the  Results 

The  solution  of  the  governing  equations  gives  the  displacements,  velocities,  and 
accelerations  of  the  system’s  various  masses.  These  results  can  be  analyzed  for  possible 
design  implications.  Accuracy  of  a  model  is  determined  by  how  well  the  solution  of  the 
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governing  equations  of  motion  predict  the  observed  behaviors  of  the  system  presented  in 
experimental  results. 

1.4  Objectives 

The  objectives  of  this  study  were  as  follows: 

1 .  to  model  the  optical  platform  as  a  rigid  plate  supported  by  elastic  springs  containing 
discrete  masses 

2.  to  determine  an  exact  differential  of  motions  representing  the  kinematic  model 

3.  to  determine  its  mechanical  response  for  free  vibration,  impact  analysis,  and  forced 
vibration. 

4.  to  incorporate  MATHEMATICA  as  the  computation  tool  for  analysis  of  the  mechanical 
response-based  method  to  analytically  solve  for  the  equations  of  motion 
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CHAPTER  2 
MODELING 

2.1  Introduction 

This  chapter  introduces  a  brief  overview  of  methods  available  to  develop  the  equation  of 
motions  for  rigid  point  supported  plates  in  order  to  provide  a  framework  of  vibration;  that  is,  to 
represent  the  physical  problem  mathematically.  Differential  equations  governing  the  forced 
vibration  of  a  rigid  point  supported  plate  were  then  developed  using  Lagrange’s  equation.  By 
defining  the  plate’s  motion  in  terms  of  generalized  displacement  parameters,  one  captured  all  six 
degrees  of  freedom  in  terms  of  three  translations  and  three  rotations. 

Methods  to  represent  or  model  physical  problems  all  have  their  roots  in  dynamic 
equilibrium  or  Newton’s  2nd  Law.  For  rigid  bodies,  two  vectors  equations  given  as 

y  F  =  m  aG 

YjMG=Ip(Xp  (2.1) 

where  p  takes  on  indices  of  x,  y,  or  z.  Here  aG  is  the  acceleration  at  the  center  of  mass,  m  is  the 

mass,  Ip  is  the  mass  moment  of  inertia  about  the  plh  axis,  and  a  the  angular  acceleration  about  the 
pth  axis.  As  the  system  became  more  complicated,  establishing  the  vector  required  relationships 
in  Newton’s  approach  became  increasingly  difficult  and  required  a  more  concise  approach. 

As  an  alternative  to  this  vector-based  formulation,  the  method  of  principle  virtual  work 
was  a  scalar-based  approach  providing  a  shift  in  formulation  from  a  Newtonian  to  a  Lagrangian 
approach.  Problems  in  dynamics  and  vibration  required  introducing  d’Alembert  principle,  since 
the  principle  of  virtual  work  was  essentially  a  variation  principle  for  problems  in  statics 
coordinates  used  to  define  the  degrees  of  freedom  for  the  system.  Stated  simply 
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The  total  work  performed  by  the  effective  forces  through  infinitesimal  virtual 
displacements  compatible  with  the  system  constrains  is  zero. 
and  mathematically  is  represented  as 


YJ(F,-pI)-Sr,  =0 


(2.2) 


In  equation  2.2,  8rt  is  a  virtual  displacement  compatible  with  the  constraints  on  the  system. 
When  deriving  equations  of  motion,  independent  coordinates  were  preferred  over  system 
position  variables,  which  lead  to  the  definition  of  generalized  coordinates  and  Lagrange’s 
Equations  of  Motion.  A  complete  derivation  of  Lagrange’s  equations  can  be  found  in  [91.  Here 
it  is  stated  as 


d_ 

dt 


fdT^ 


'  })cl,  j 


dT  dR  dV  „ 

- + — +  —  =  Qi 

dcji  dqj  dqt 


(2.3) 


where  T  is  the  system’s  kinetic  energy,  V,  the  system’s  potential  energy,  and  R  is  defined  as 
Rayleigh’s  dissipation  function  and  represents  energy  loss  due  to  damping.  Also  q;’s  are 
generalized  displacements  and  Q;’s  are  the  corresponding  generalized  forces.  The  term 
generalized  is  used  to  described  translations  or  rotations  as  well  as  forces  or  moments. 


2.2  Kinematic  Model  Development 

Lagrange’s  equation  of  motion  was  the  preferred  approach  to  develop  the  differential 
equations  of  motion  for  the  point  supported  plate.  To  do  so  required  a  definition  of  the 
generalized  displacements,  and  provided  the  starting  point  of  kinematic  analysis. 

Ligure  2.1  shows  the  experimental  setup  of  the  optical  platform  and  all  of  its  required 
instrumentation.  Two  inertial  actuators  were  mounted  to  provide  forced  motion  in  the  vertical 
direction  as  well  as  to  create  angular  rotations  about  the  three  cartesian  axes.  As  a  consequence 
of  the  induced  motions,  displacements  in  the  x  and  z  directions  were  produced  and  accounted  for 
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in  the  kinematic  model.  Other  elements  included  mounted  brackets  for  optical  sensors,  laser 
beam,  and  sensors  themselves.  Each  element  added  mass  and  hence  affected  the  mechanical 
response  of  the  optical  platform  system.  Vibration  isolators  provided  a  mounting  base  for  each 
of  elastic  four  springs. 


Mount  for  Sensor 


Inertial  Actuator 


Mountfor  Sensor 


Isolatorand  Spring 


Optical  Platform  Inertial  Actuator 


Figure  2.1:  Experimental  Setup 

As  such,  the  system  was  modeled  as  an  elastically  point  supported  plate  containing  discrete 
masses.  The  platform  was  taken  to  be  a  rigid  body,  which  simulated  the  construct  of  the  actual 
system  in  operation. 

The  developed  kinematic  model  consisted  of  superimposed  displacement  in  the  x,  y,  and 
z  directions.  Three  two-dimensional  representations  were  used  to  develop  the  total  displacement 
expressions  u{y,z},  v{x,z}and  w{x,y}  for  the  system  consisting  of  a  rigid  platform  and  any 
number  of  discrete  masses,  depicted  below  as  the  klh  masses.  The  cartesian  coordinate  system 
had  been  located  at  the  geometric  center  of  the  plate,  its  centroid,  and  provided  the  means 
necessary  to  establish  a  relative  position  vector  expression  to  locate  movement  in  terms  of 
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displacements  given  by  {uo,  Vo,  wo}  and  angular  rotations  given  by  { 0X,  0y,  0Z}.  The  procedure 
used  to  develop  the  displacement  expression  is  presented  in  equations  (2.4)  through  (2.9).  First 
consider  figure  2.2,  which  shows  kinematics  relationships  for  the  x-y  plane.  All  terms  were 
considered  positive.  The  centroid  displaced  by  Uo  and  Vo,  which  represented  the  translation  of 
the  center  of  gravity.  In  this  plane,  0Z  represented  a  positive  rotation,  and  elements  located  in  the 
positive  x-y  quadrant  experienced  displacements  given  by 

rk/G  ={j c-  xCos  (0,)  -  y  Sin  (0z)}i  +  {x Sin  (0:)  -  (y  -  y  Cos  (0.))}  ]  (2.4) 

Combining  translations  and  rotations  provided 

rk  =  {m0  +  x  -  x Cos  (6, )  -  y  Sin  (6. ) } i  +  {v0  +  x Sin  (&z)  —  (y  —  y  Cos  {6Z )) }  j  @-5) 


k*  mass 


Figure  2.2:  Free  Body  Diagram  in  X-Y  Plane 

The  amplitude  of  motion  experienced  by  the  plate  and  its  discrete  elements  were  on  the  order  of 
micrometers  and  micro-radians.  Using  a  small  angle  assumption  resulted  in 


h  ={uo-y0z}i  +{vo+x(0z)}j 


(2.6) 
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Referring  to  figure  2.3  and  making  similar  assumptions  of  displacement  and  rotation  amplitudes, 
the  position  vector  in  the  y-z  plane  was  computed  as 

rk  ={vo-z0x}]  +  {wo  +  y  6X } k  (2.7) 

The  final  contribution  was  determined  by  referring  to  figure  2.4,  which  represented  motion  in  the 
x-z  plane.  This  became 

rk  =  {u0  +  x +  z  9y }  I  +  {w0  -  x  9y  +  z }  k  ^8) 

The  total  position  vector  now  was  determined  by  summing  equations  (2.5)  to  (2.8)  producing 

rk  =u{y,z,t}i  +v{x,z,t}  j  +  w{x,y,t}k 

where 

u(y,z,t)  =  u0(t)+ z9y(t)~  y9z  (t) 
v(x,z,t )  =  v0(t)~  z9x(t)  +  x0z(t) 
w(x,  y,  t )  =  w0  (t)  -  x9y  it)  +  y9x  (t) 

Care  was  taken  to  not  account  for  contributions  more  than  once. 


(2.9) 

(2.10) 


Figure  2.3:  Free  Body  Diagram  in  Y-Z  Plane 
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Equation  (2.10)  is  the  displacement  equation  required  for  Lagrange’s  equation. 


W0 


'  r 


Figure  2.4:  Free  Body  Diagram  in  X-Z  Plane 

2.3  Dynamic  Model  Development 

Figure  2.5  shows  a  planar  view  of  dynamic  model  that  represented  the  optical  platform. 
The  actual  system  had  two  mounted  actuators  oriented  45°  relative  to  the  horizontal.  Each 
inertial  actuator  had  been  modeled  as  an  equivalent  mass-spring-damper  system  with  stiffness  ksj 
and  damping  coefficient  csj .  Here  mtopwas  the  portion  of  the  actuator’s  mass  in  motion  and 
mbase  was  the  remaining  mass  for  the  actuator.  Discrete  masses  were  represented  by  the  kth  mass. 
Experimental  testing  of  the  system’s  support  springs  revealed  that,  although  geometrically 
similar,  measured  spring  stiffnesses  could  vary  by  as  much  as  30%.  To  account  for  this 
variation,  spring  stiffnesses  and  damping  coefficients  were  all  labeled  differently  as  kj  and  c, . 
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Figure  2.5:  Side  view  of  a  Platform  with  Discrete  Mass  and  Actuator 
Expression  for  the  kinetic  energy,  potential  energy  and  Rayleigh  Dissipation  function  for 


each  inertial  actuator  were  given  by 


=  J  msX  («A  +  ^  )  +  J  ms2  +  VV,22  ) 


(2.11) 


^ac, valors  =  1  K“,l  “  U  (  V  A  .  z  U  »' '  +  (VU  “  W  ?,! »  Z,1  1)  "  +  M  ?,1 >  Z,1  1)1  '  y  ^ 

+  -ks2\-u{yS2’Zs2}2  +(vS2-v{ys2,zs2})2  +(ws2-w{ys2,zs2})]2 

z  | 

=-CJi[(M,i-M{y,1,Zsi})2  +(vsl  -vtyri.Z,!))2  +vv{ysl,zil})2] 


2 


+  -csi{u{ys2iZs2)2  +  (vi2  ~v{ys2,zs2})2  +  (ws2  -w{ys2,zs2})2}  (2.i3) 


and  expressions  for  the  kinetic  energy,  potential  energy,  and  Rayleigh’s  Dissipation  function  for 
the  plate  were  given  as 

T  +  D*  +Wi2)  +  ^\\\(Vi2  +E2  +w2)dm 

^  i= 1  ^ 

R  =  ^c,(m,2+v;2+w,2) 

^  1=1 

V  =^^,(m,2  +v;2  +w;2) 

^  i=l 


(2.14) 
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Substitution  of  equations  (2.12)  through  (2.15)  into  (2.3)  yielded  the  coupled  equations  of  motion 
represented  in  matrix  form  as 


1 

Co 

4^ 

X 

4^ 

[0](4.v6) 

| {<ls  }(4xl)  |  + 

r  rc  l 

Lv  1(4x4) 

IT/as  ]  (4y6) 

SO 

o 

_ 1 

K](  6x6). 

sp  ](6jc4) 

\-C p  ](6y6) 

+ 

[^j](4y4) 

](4.v6) 

[ksp  ](6.v4) 

[k p  ](6.v6) 

fe.i-  }(4.vl)  1 

P  J(6y1)  J 

}(4.vl)  | 

p  j"(6y!)  J 


[{&}<  4x1)1 

|  {0W)  J 


(2.15) 


The  explicit  equations  were  quite  lengthy  and  can  be  referenced  in  Appendix  A  - 1  and  appendix 
A-2,  which  provide  sub-matrices  used  in  eqn.(2-15).  The  above  subscripts  s  and  p  in  eqn.  (2.16) 
correspond  to  the  shaker  and  plate,  respectively.  Also,  the  vectors  of  generalized  coordinates  are 
defined  as 


kl= 


s2 


W, 


s  2. 


U= 


w0 

8, 

e. 


(2.16) 


This  set  of  equations  is  coupled  both  statically  and  dynamically,  representing  damped  forced 
vibration  analysis.  Initial  condition  must  be  appended  to  eqn  (2.15),  which  consisted  of  initial 
displacement  and  initial  velocities. 

Equation  (2.15)  was  specialized  based  upon  the  type  of  vibration  analysis  investigated. 
Initial  conditions  can  also  be  presented.  For  free  damped  vibration  analysis,  eqn.  (2-15)  reduced 
to 


\-mp  ](6y6)  p  ]( 6x6)  \flp  1(6*1)  ^  p  ](6y6)  \(,x\)  {^J(6xl) 


(2.17) 
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with  initial  conditions  of 


M0(0) 

o' 

v0(0) 

0 

W’c,  (0) 

fed- 

0 

*,(0) 

0 

*<(0) 

0 

A-(0). 

0 

V.  .7 

(2.18) 


Although  initial  velocity  was  possible,  only  displacement  conditions  were  presented  to  reflect  the 
reality  of  inducing  such  initial  conditions  in  a  laboratory  setting.  Equation  (2.17)  also  governed 
impact  analysis.  However,  the  accompanying  set  of  initial  conditions  were 


fed- 


0 

0 

0 

V0(0+) 

0 

fed- 

0 

0 

A(o+) 

0 

0 

0 

A-(0+). 

These  initial  conditions  were  computed  from  changes  in  both  linear  and  angular  momentum  of 
the  plate  just  before  impact  and  immediately  after  impact,  given  by 


o+ 

|  Fdt  =  mv(0+ )  -  mv(0~ ) 


0“ 


o+ 


0“ 

0+ 


0“ 


(2.19) 


The  impact  location  did  have  an  effect  on  the  initial  conditions.  A  plate  impacted  in  its  center 
resulted  only  in  changes  in  linear  momentum  and  not  angular  momentum. 
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Finally,  steady-state  forced  vibration  analysis  for  the  system  was  considered.  No  initial 
conditions  were  required. 
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CHAPTER  3 
MODAL  ANALYSIS 

3.1  Modal  Analysis 

In  this  chapter  an  overview  of  the  solution  strategy,  modal  analysis,  is  presented  and  used 
to  solve  the  mechanical  response  in  light  of  the  complicated  structure  of  the  governing  equations. 
The  equations  of  dynamic  equilibrium  were  coupled  both  statically  and  dynamically  and,  with 
the  presence  of  damping,  did  not  render  a  decoupled  system  unless  certain  conditions  existed 
between  the  stiffness,  damping  and  mass  matrices. 

There  were  many  ways  to  decouple  eqn.  (2.15)  based  upon  the  existing  type  coupling 
construct.  Consider  the  two-degree  of  freedom,  coupled  system,  given  by 


or  in  matrix  form  as 


(3.1) 


[m]{x}  +  [c]{i}  +  [fc]{4  =  {L}  (3.2) 

Now  simplify  the  equations  to  consider  an  undamped  problem  with  mi2  =  rmi  =  0,  a  system 
coupled  by  virtue  of  the  stiffness  elements  only.  Since  the  mass  matrix  is  now  diagonal,  a 
coordinate  transformation,  defined  as 


i 

{q}  =  [m]2{x} 

produces 


(3.3) 


i  _i 

[m][mr-{q}+[k][m]^{q}  =  {F} 


(3.4) 


_i 

Pre-multiplying  eqn. (3. 3)  by  [ m ]  2  yields 


[I]{q}+[k]{q}  =  {F} 


(3.5) 
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where 


i  i 

[/]  =  [m]  2[m][m ]  2 


i  i 

[k]  =  [m]  2[k][m]  2 


_i 

{F}  =  \m]2{F} 


(3.6) 


In  eqn.  (3.5),  | k  |  is  the  mass  normalized  stiffness  matrix  and  [I]  is  the  identity  matrix.  The 
modal  transformation  matrix  is  accomplished  by  determining  the  matrix  of  eigenvectors 
corresponding  to  eqn.  (3.5),  identified  here  as  [P],  defining  a  modal  coordinate  system  {r} 
through  the  transformation  of  {q}  =  [P]{r},  and  incorporating  this  transformation  into  eqn.  (3.4). 
Doing  so  produces 

mm{r}+mm{r}  =  {F}  (37 

T 

and  pre-multiplying  by  fP]  yields  a  decoupled  system  of  equations  of 


{f}  +  [k]{r}  =  {F} 

where 

[k]  =  [P]r[k][P] 
{F}  =  [Pf{F} 


(3.8) 


(3.9) 


Eqn  (3.8)  was  solved  directly.  The  solution  to  the  required  problem  was  obtained  after  the 
solution  in  modal  space  was  transformed  back  to  the  original  coordinate  system  using  the  modal 
transformation  given  by 

W  =  [5]W  (3.10) 

where 

i 

[S]  =  [m]2[P]  n  in 


Now  consider  the  un-damped  problem,  but  one  that  is  coupled  both  statically  and 
dynamically;  that  is,  one  with  fully  populated  mass  and  stiffness  matrices.  Using  the 
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transformation  defined  through  eqn.  (3.3)  is  not  the  most  efficient  method  for  this  type  of 
problem.  Rather,  a  Cholesky  Decomposition  is  the  preferred  strategy. 

The  Cholesky  Decomposition,  L,  produces  a  lower  triangular  matrix  that  satisfies  the 

T 

matrix  product  of  [m]  =  |  L||  L|  .  Most  of  the  context  presented  will  hold  for  this  coordinate 
transformation,  defined  as 


which  produces 


{q}  =  [LT]{x} 


[m][LTrl{q}  +  [k][LT]l{q}  =  {F} 


Pre-multiplying  eqn. (3. 7)  by  \L\  1  yields 

U\{q}  +  [k]{q}  =  {F} 


where 

[I]=[L]-l[m][LT]-1 
[k]  =  [LYl[k][LT  ]-' 
{F}  =  [L]  ‘{F} 


(3.12) 

(3.13) 

(3.14) 

(3.15) 


The  eigenvector  matrix  [P]  can  now  be  determined,  and  the  modal  transformation  matrix  [S] 
required  is  defined  through 

M  =  [S]M  (3.16) 

where 

[5]  =  [Lr]-,[P]  (3.17) 

[sr1=[pr][Lr] 

Damping  complicated  the  process  in  several  ways.  First,  it  was  challenging  to  measure,  and 
secondly,  with  damping  present,  there  were  no  assurances  that  the  matrix  of  eigenvectors  would 
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diagonalize  the  damping  matrix,  thus  decoupling  the  system.  Caugley  and  O’Kelly  [13]  have 
shown  exactly  that. 


[k][mY1[c]  =  [c][m]l[k] 

[m][kY1[c]  =  [c][kY1[m]  (3-18) 

[m][c]-l[k\  =  [k][c]-\m] 

Whenever  eqn.  (3.18)  does  not  hold,  then  a  proportional  damping  model  is  recommended  in 
which  the  damping  matrix  [c]  can  be  expressed  as  a  linear  combination  of  the  mass  and  stiffness 
matrix  gven  by 


[c]  =  a[M]  +  J3[k]  (3.19) 

Constants  a  and  [3  were  selected  to  produce  a  desired  damping  ratio,  based  upon  experimental 
results  or  design  considerations. 

Since  the  damping  is  mode  dependent,  the  damping  ratio  is  given  as 


c, 


a  ,  Poii 
2(0,  2 


(3.20) 
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CHAPTER  4 
FREE  VIBRATION 

4.1  Introduction 

This  chapter  presents  the  free  vibration  analysis  of  the  system.  Equation  (2.17)  and 
(2.18)  were  solved  based  upon  the  discussion  in  chapter  3.  Required  in  the  analysis  were  all  of 
the  system  properties.  The  mass  or  the  platform  and  discrete  elements  were  weighted  and 
recorded.  Data  from  the  manufacturer  of  the  support  springs  indicated  the  same  spring  stiffness. 
However,  experimental  results  for  the  spring  stiffnesses  indicated  stiffness  coefficients  that 
varied  as  much  as  30  percent.  Free  vibration  case  studies  were  investigated,  which  vary  spring 
stiffness,  discrete  mass  amount  and  location,  and  in  the  case  of  damping,  damping  ratios.  Target 
displacement  plots  were  plots  that  map  the  plate’s  motion  onto  a  plane  located  a  set  distance 
away  from  the  center  of  the  plate,  and  presented  a  visual  viewpoint  of  the  platform’s  local 
motion.  All  results  computed  herein  were  performed  using  MATHEMATICA®. 

Figure  4.1  shows  the  topology  of  stiffness  and  masses  used  in  the  vibration  studies.  The 
coordinate  system  was  placed  at  the  centroid.  The  plate  had  geometry  of  36  in  x  36  in  x  2  in. 
with  a  density  of  0.016  lb/in  .  The  springs  were  mounted  at  fixed  coordinates  as  listed  in  Table 


4.1. 
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Figure  4.1:  Mass,  Stiffness,  and  Damping  Coefficient  Distribution 

Parametric  evaluation  of  the  effect  of  spring  stiffness,  discrete  mass  composition  and  discrete 
mass  location  values  were  established  in  Tables  4.1  through  4.4.  Table  4.1  shows  the  fixed 
location  of  the  support  springs  while  Table  4.2  presents  the  spring  stiffnesses  used.  Table  4.3 
presents  the  discrete  masses  used  in  the  study,  and  Table  4.4  shows  how  the  masses  were 
distributed  on  the  plate.  A  systematic  notation  to  identify  what  type  of  vibration  was  considered, 
the  value  of  stiffnesses  used,  the  location  of  the  discrete  masses,  and  the  amount  of  each  discrete 
mass  is  given  as 

z  -  £/  -  -  U 

Vibration  Type  Stiffness  Case  Mass  Case  Mass  Location 

where  V  identifies  the  type  of  vibration  problem  as  FR  for  free  vibration,  IM  for  impact  analysis 
and  FF  for  forced  vibration.  For  example,  IM-K1-M2-L3  identifies  an  impact  analysis  using 
equal  spring  stiffnesses  with  equal  masses  distributed  randomly  over  the  plate’s  area.  In  the  case 
of  no  added  mass,  the  mass  location  index  was  omitted  and  would  append  as  V-K,-Mj. 
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Table  4.1  Spring  Coordinate  Location 

Spring 

Number 

Cartesian  Coordinates 

x  (in) 

y  (in) 

z  (in) 

1 

17.25 

-1 

17.25 

2 

-17.25 

-1 

17.25 

3 

-17.25 

-1 

-17.25 

4 

17.25 

-1 

-17.25 

Table  4.2  Spring  Stiffnesses 

Spring 

Stiffness  (lb/in) 

Case  Study 

1 

2 

3 

4 

5 

ki 

115 

115 

115 

110 

115 

k2 

115 

110 

115 

110 

110 

k3 

115 

110 

110 

115 

120 

k4 

115 

115 

110 

115 

100 

Table  4.3  Discrete  Masses 

Mass  (slugs) 

Case  Study 

1 

2 

3 

4 

5 

6 

mi 

0.0 

0.155 

0.155 

0.155 

0.078 

0.155 

rm 

0.0 

0.155 

0.078 

0.155 

0.155 

0.078 

m3 

0.0 

0.155 

0.078 

0.078 

0.078 

0.0 

m4 

0.0 

0.155 

0.155 

0.078 

0.155 

0.233 

Table  4.4  Discrete  Mass  Location 

Location 

Case  Study 

(Xl,Zi) 

(x2,z2) 

(x3,z3) 

(x4,z4) 

li 

(9,9) 

(12,12) 

(12,12) 

(9,9) 

h 

(-9,9) 

(6,6) 

(12,6) 

(-12,12) 

13 

(-9,-9) 

(12,-6) 

(6,6) 

(-15,-5) 

u 

(9,-9) 

(6,-12) 

(6,12) 

(6,-6) 
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4.2  Free  Vibration  Analysis 

Free  vibration  analysis  was  governed  by  eqns.  (2.17)  and  (2.18).  Results  consisted  of 
natural  frequencies,  displacement  and  rotation  results,  and  target  plots  for  various  types  of  initial 
conditions.  Initial  conditions  consisted  of  non-zero  generalized  displacements  and  zero 
generalized  velocities. 

Figure  4.2  shows  results  for  displacements  and  rotations  for  FR-KrM|  with  an  initial 
displacement  of  0.01  ft  in  the  x-direction  only.  Natural  frequencies  of  this  case  were 
coi  =  C02=  5.464  Hz,  003=  5.478  Hz,  0)4=  9.094  Hz,  0)5=  C06  =  9.104  Hz.  Modes  1  and  2 
correspond  to  the  in-plane  modes,  mode  3  is  the  transverse  mode,  mode  4  corresponded  to  the 
yaw  mode,  and  modes  5  and  6  corresponded  to  the  rotational  modes.  The  non-zero  displacement 
was  u0(t)  and  the  only  non-zero  rotation  was  0Z{ t) .  These  expressions  were 

Uq (t)  =  2.967X  10  5  Sin[1.5708  -  57.203 1]  +  9.970x  10  3  Sin[1.5708  +  34.3343 1]  (4.1) 

6Z  (t)  =  6.270 x  10  4  Sin[  1 .5708  -  57.203 1]  -  6.270 x  10’4  Sin[  1 .5708  +  34.3343 1] 


Both  motions  were  in-phase  with  each  other. 
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Figure  4.2  Displacement  plots  for  free  vibration  equal  stiffnesses,  no  added  masses, 
with  initial  condition  of  uo(0)  =  0.01  ft. 

Changing  the  stiffness  case  to  case  5,  or  evaluating  FR-K5-M1  for  the  same  initial  condition, 
highlighted  the  effect  of  unequal  spring  stiffness  and  static  coupling  on  the  motion  of  the 
platform.  All  displacements  and  angular  rotations  grew  harmonically  with  the  exception  of 
Uo(t).  This  induced  motion  exhibited  a  beating  phenomena  characterized  by  a  slow  variation  in 
amplitude  with  rapid  oscillations. 

Clearly  based  upon  spring  stiffnesses,  two  types  of  motions  exist.  The  first  is  identified 
as  direct  motion  consistant  with  initial  conditions  and  the  second  is  indused  motion  based  upon 
system  properties.  Another  noticeable  result  was  that  the  natural  frequencies  separated  to  coi  = 
5.402  Hz,  co2  =  5.404  Hz,  co3=  5.418  Hz,  co4=  8.804  Hz,  ro5=  8.996  Hz,  co6=  9.202  Hz. 
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Figure  4.3  Displacement  plots  for  free  vibration  unequal  stiffnesses,  no  added 
with  initial  condition  of  uq(0)  =  0.01  ft. 
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Expressions  for  the  displacements  highlighted  the  coupling  effects  and  were  given  as 
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u0(t)  =  1.702  xlO'5  Sin[1.5708 -55.3206  t]  + 1.281x10  3  Sin[l. 5708  -  34.084  t] 

+  6.288  x  10'3  Sin[1.5708  +  33.9414  t]  +  2.399x  10  3Sin[1.5708  +  33.9563  t] 
+  3.120x  10'8Sin[1.5708  +  56.52041 1]  +  1.300x  10  5Sin[1.5708  +  57.8158  t] 

Vg ( t )  =  -  6.525  x  10  6  Sinfl .5708  -  55.3206  t]  +  3.342  x  10'3  Sinfl .5708  -  34.084  t] 
-2.410X10'3  Sin[1.5708  +  33.9414  t]-9.199x  10  4Sin[1.5708+ 33.9563  t] 
-1.196x10  sSin[1.5708+ 56.52041  t]-4.984x  10'6Sin[1.5708  + 57.8158  t] 

w0(0  =-1.692xl0'5  Sin[1.5708- 55.3206 1]  - 4.158x10  3  Sin[1.5708  +  33.9414  t] 
+  4. 162  x  10‘3  Sin[1.5708  +  33.9563  t]  -  2.639  x  10  7  Sin[1.5708  +  56.5204  t] 
+ 1.307 xl0  5Sin[l. 5708  +  57.8158  t] 

ex(t)  =3.363 x  10  4  Sin[1.5708  -  55.3206  t]  - 2.940 x  10‘4  Sinfl .5708  +  33.9414  t] 
+  2.508  x  10'4  Sin[l .5708  +  33.9563  t]  +  3.440 x  10  8Sin[l .5708  +  56.204  t] 

-  2. 93 lx  10  4 Sin[  1.5708  +  57.8158  t] 

(0  =  8.299  x  10‘8  Sin[l  .5708  -  55.3206  t]  +  4.997  x  I0'5  Sinfl  .5708  +  33.9414  t] 

-  5.009x  10  5  Sinfl.5708  +  33.9563  t]  -  1.441x  10  5Sin[1.5708  +  56.204  t] 

+  6.232  xl0'6Sin[1.5708  +  57.8158  t] 

6z(t)  =  3.349 XlO’4  Sin[1.5708- 55.3206  t]- 4.737x10  4  Sin[l. 5708  +  33.9414  t] 
- 1 .557  x  10  4  Sin[l .5708  +  33.9563  tl  +  6.573  x  10‘7Sin[1.5708  +  56.5204  t] 
+  2.938  x  10"4Sin[1.5708  +  57.8158  t] 


Figure  4.4  shows  the  results  for  initial  conditions  of  Uo(0)  =  0.01  ft,  Vo(0)  =  -0.01  ft  and  Wo(0)  = 
0.01  ft,  with  all  other  initial  conditions  equal  to  zero,  since  it  is  more  difficult  to  prescribe 
angular  position  than  displacements.  Here  notice  that  Uo  and  Wo  move  in-phase  while  0y  and  0Z 
move  in-phase. 
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Figure  4.4  Displacement  plots  for  free  vibration  unequal  stiffnesses,  no  added  masses 
with  initial  condition  of  u0(0)  =  0.01  ft,  v0(0)  =  -0.01  ft,  wo(0)  =  0.01  ft 

Below  are  target  plots,  which  represent  a  two-dimensional  projection  of  the  motion  of  the 

platform  and  gave  an  indication  of  how  the  motions  increased  with  distance  and  provided  a 

visual  perspective  of  the  platform’s  motion.  Here  the  distance  was  14  ft  from  the  center  of  the 

plate.  Figure  4.5a  corresponded  to  FR-Ki-Mi  and  figure  4.5b  corresponded  to  FR-K5-M1. 
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Figure  4.5  Target  Plots  for  equal  and  unequal  spring  stiffness  and  no  added  discrete  mass. 


Adding  mass  reduced  the  natural  frequency  of  the  system,  as  well  as  the  displacement 
amplitudes.  Increasing  the  mass  at  points  not  co-located  with  the  centroid  affected  rotational 
modes  of  vibration.  Take,  for  example,  the  case  identified  as  FR-KrM2-L| ,  which  added  equal 
discrete  mass  of  0.155  slugs  at  the  center  of  each  quadrant.  The  natural  frequencies  were  coi  =  cco 
=  5.114Hz,  C03  =  5.147  Hz,  a>4  =  005=  9.104  Hz,  and  cc>6=  8.672  Hz.  Modes  1  and  2 
corresponded  to  the  in-plane  modes,  mode  3  was  the  transverse  mode,  mode  4  and  5  were  the 
rotational  modes,  and  mode  6  was  the  yaw  mode.  Here  there  was  a  shift  in  higher  rotational 
modes.  An  interesting  scenario  was  the  case  of  unequal  stiffness,  with  added  unequal  discrete 
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masses  not  arranged  in  a  symmetric  pattern  about  the  origin,  as  depicted  by  FR-K5-M6-L4.  This 
resulted  in  natural  frequencies  that  were  all  unique  given  by:  ©1  =  5.047  Hz,  ©2  =  5.055Hz,  ©3  = 
5.0885  Hz,  ©4=  8.324  Hz,  ©5=  8.498  Hz,  and  ©6=  8.557  Hz.  Displacement  plots  are  shown  in 
figure  4.5,  and  its  target  plot  is  presented  in  4.6. 
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Figure  4.6  Displacement  plots  for  free  vibration  unequal  stiffnesses,  added  unequal  discrete 
masses  with  initial  condition  of  vo(0)  =  -0.01ft 
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Figure  4.7  Target  Plot  for  unequal  spring  stiffness  and  added  discrete  mass. 

4.3  Free  Vibration.  Damped  Analysis 

In  this  section,  damping  was  included  in  the  vibration  analysis  and  represented  using  a 
proportional  damping  model  as  presented  in  eqn.  (3.20).  Initial  conditions  selected  were  similar 
to  those  used  previously.  Several  combinations  of  constants  a  and  (3  were  computed  in  order  to 
yield  realistic  damping  ratios.  Results  presented  in  this  section  followed  a  pattern  similar  to  the 
results  presented  in  the  previous  section. 

Table  4.5  compares  values  of  natural  frequency  con,  damped  natural  frequency  cod,  and 
damping  ratio  C,  based  upon  proportional  damping  constants  a  and  [3.  Damping  ratios  of  the 
order  of  1  to  2  %  are  realistic  and  provide  a  framework  for  selecting  proportional  constants.  For 


convenience  these  constants  were  taken  to  be  0.0007812  for  damped  free  vibration,  although  any 
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set  can  be  chosen  based  upon  performance  and  design  considerations.  These  values  provided  a 
maximum  damping  ratio  of  2.23%. 


Table  4.5  Natural  Frequency,  damped  natural  frequency  and  damping  ratio  as  function  of 

proportion  constants  a  and  (3 


Mode 

a  =  0.0125 
(3  =  0.0125 

a  =  0.00625 
(3  =  0.00625 

a  =  0.003125 
(3  =  0.003125 

C0n 

cod 

c 

C0n 

cod 

c 

C0n 

cod 

c 

1 

5.464 

4.842 

0.215 

5.464 

5.163 

0.107 

5.464 

5.315 

0.0537 

2 

5.464 

4.842 

0.215 

5.464 

5.163 

0.107 

5.464 

5.315 

0.0537 

3 

5.478 

4.853 

0.215 

5.478 

5.176 

0.108 

5.478 

5.329 

0.0538 

4 

9.094 

7.291 

0.357 

9.094 

8.242 

0.179 

9.094 

8.678 

0.0893 

5 

9.104 

7.291 

0.357 

9.104 

8.250 

0.179 

9.104 

8.687 

0.0894 

6 

9.104 

7.297 

0.357 

9.104 

8.250 

0.179 

9.104 

8.687 

0.0894 

Mode 

a  =  0.001562 
(3  =  0.001562 

a  =  0.0007812 
(3  =  0.0007812 

a  =  0.000390 
(3  =  0.000390 

C0n 

cod 

c 

C0n 

cod 

c 

C0n 

cod 

c 

1 

5.464 

5.390 

0.0268 

5.464 

5.427 

0.0134 

5.464 

5.446 

0.0067 

2 

5.464 

5.390 

0.0268 

5.464 

5.427 

0.0134 

5.464 

5.446 

0.0067 

3 

5.478 

5.404 

0.0269 

5.478 

5.441 

0.0134 

5.478 

5.460 

0.0067 

4 

9.094 

8.888 

0.0446 

9.094 

8.992 

0.0223 

9.094 

9.043 

0.0111 

5 

9.104 

8.898 

0.0447 

9.104 

9.001 

0.0223 

9.104 

9.053 

0.0111 

6 

9.104 

8.898 

0.0447 

9.104 

9.001 

0.0223 

9.104 

9.053 

0.0111 

The  presence  of  damping  in  the  vibration  model  resulted  in  a  characteristic  equation,  which 
contained  complex  roots.  The  displacement  expressions  presented  in  eqn.  (4.4)  reflected  this  fact 
and  were  used  to  generate  figure  4.8  for  the  case  of  FR-Ki-Mj.  Here  the  plate  was  given  an 
initial  displacement  in  the  x  direction  only.  Damped  natural  frequencies  were  given  as  C0di  =  C0d2 


=  5.437  Hz,  C0d3  =  5.442  Hz,  C0d3  =  8.992  Hz,  C0d5  =  C0d6  =8.708  Hz.  Figure  4.9  is  the  target  plot 
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and  shows  that  the  density  of  the  plot  trace  became  more  concentrated  near  its  origin  as  time 
increased,  which  reflected  the  effect  of  damping. 


u0(t)  =2.967  xlO'5  e  1  279t  Sin[l. 5934  -  56.560  t] 

+  9.970x10  3  e  04608tSin[1.5578  +  34. 103 It]  (4.4) 

0z(t)  =  6.270xl0~4  e  L279t  Sin[  1. 5934  - 56.560 1] 

-6.270x10  4  e  04608t  Sin[1.5578  +  34.1031 1] 


Figure  4.8  Displacement  plots  for  free  vibration  unequal  stiffnesses,  no  added  masses  with  initial 

condition  of  uq(0)  =  0.01ft. 
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Figure  4.9  Target  plot  for  equal  stiffness,  no  added  discrete  mass  damped  vibration. 

If  the  spring  stiffness  case  changed  to  case  5  or  FR-K5-M1.  the  solution  for  the  six  generalized 
displacements  were  given  in  eqns.  (4.5)  and  (4.6).  The  corresponding  damped  natural 
frequencies  were  all  distinct,  with  values  of  codi  =  5.366  Hz,  C0d2  =  5.368  Hz,  C0d3  =  5.382  Hz,  cc>d4 
=  8.708  Hz,  cod5  =  8.895  Hz,  cod6  =  9.907  Hz. 


u0(t)  =  1.702xl0"5  e"1196t  Sin[1.5932-54.7195t]  +  1.281xl0"3  e"°'453t  Sin[  1 .5842- 33.82 lit] 

+  6.288xl0"3  e"°'4503t  Sin[1.5574+33.7155t]  +  2.399xl0'3  e  04507t  Sin[  1 .5574+  33.7301 1] 
+  3.120xl0"8e_1'2482t  S  in[l. 5484+55. 8928t]  +  1.300xl0"5  ea4507tSin[1.5479+57.159t] 

v0(f)  =  -6.525xl0'6  e'U96t  Sin[1.5932-54.7195t]  +  3.342xl0'3  e  a453t  Sin[  1 .5842- 33.82 1 1 1] 
-2.410xl0"3  e"°'4503t  Sin[l .5574+  33.7 155 1] -9. 199x1 0"4e"°'4507t  Sin[1.5574+33.7301t] 
— 1.196xl0'8  e"L2482t  Sin[  1.5484+ 55. 8928 1]- 4. 984xlO"6e"a4507tSin[  1.5479+ 57. 159 1] 

w0(t)  =-  1.692x1 0‘5  e"1196t  Sin[1.5932-54.7195t]-4.158xl0-3  e“a453t  Sin[1.5574+  33.7155 1] 

+  4.162xl0"3  e"°'4507t  Sin[l. 5574+33. 7301t]-2.639xl0"7  e  L2482t  Sin[l. 5484+ 55.8928 1] 
+ 1 ,307x  1 03e  045(,7lSin[  1 .5479+  57. 1 59 1] 


(4.5) 
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dx(t)  =3.363xl0"4  e"U96t  Sin[l. 5932- 54.7 195t]-2.940xl0"4  e  0453t  Sin[1.5574+ 33.71551] 
+  2.508x1  O' Va4507t  Sin[1.5574+33.7301t]  +  3.440xl0VL2482tSin[1.5484+55.8928t] 

-  2.931xlO"4e"°'4507tSin[1.5479f  57.159t] 

dy(t)  =8.299xl0"8  e"U96t  Sin[1.5932-54.7195t]  +  4.997xl0"5  e  °'453t  Sin[1.5574+ 33.71551] 

-  5.009x10"  Va4507t  Sin[l. 5574+  33.730 It] -1.44 lx  10"V  24821  Sin[l. 5484+  55.89281] 
+  6.232xl0Va4507tSin[l  .5479+  57. 1591] 

e7(t)  =3.349xl0"4  e  "96t  Sin[l. 5932- 54.7 195t]-4.737xl0"4  e  0453t  Sin[1.5574f  33.71551] 
-1.557xlO"4e"°'4507t  Sin[l. 5574+ 33.73011] +  6.573xl0"7e  L2482tSin[l. 5484+ 55.89281] 
+  2.938xl0"4e"°'4507tSin[1.5479+57.159t] 

(4.6) 

Figure  4.10  shows  the  displacement  as  a  function  of  time  using  a  and  [3  =  0.0007812.  Both  uo 
and  Vo  were  more  dominant  in  respect  to  displacement  magnitudes,  since  wowas  on  the  order  of 
10"5,  while  0y,  on  the  order  of  10’6,  was  less  significant  when  comparing  the  angular  rotations. 

The  selection  of  constants  a  and  [3  was  based  upon  a  desire  to  have  a  maximum  damping 
ratio  of  approximately  2%.  Increasing  the  constants  raised  ^  to  a  maximum  of  4%  which 
dampened  out  motions  somewhat  faster  as  shown  in  figure  4.12.  There  was  a  change  in  the 
damped  frequencies  to  C0di  =  5.037  Hz,  C0d2  =  5.3039  Hz,  C0d3  =  5.052  Hz,  C0d4=  8.610  Hz,  C0d5  = 


8.795  Hz,  and  C0d6  =  8.990  Hz. 
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Figure  4.10.  Displacement  plots  for  free  vibration,  damped  unequal  stiffnesses,  no  added 
discrete  masses  with  initial  condition  of  Uo(0)  =  0.01ft  with  a  =  0.0007812,  (3  =  0.0007812 
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Figure  4.1 1  Target  Plot  for  unequal  spring  stiffness  and  no  added  discrete  mass  (damped 
vibration). 
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Figure  4.12  Displacement  plots  for  free  vibration  unequal  stiffnesses,  added  discrete 

masses  with  initial  condition  of  uq(0)  =  0.01ft  with  a  =  0.001562,  (3  =  0.001562 
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CHAPTER  5 
IMPACT  ANALYSIS 

5.1  Introduction 

Before  moving  onto  forced  vibration,  some  results  are  presented  for  impact  analysis. 
Impact  analysis  had  the  same  governing  differential  equations  as  free  vibration,  but  here  velocity 
and  angular  velocity  initial  conditions  were  prescribed  instead  of  displacements  and  angular 
rotations.  These  conditions  corresponded  to  a  change  in  both  linear  and  angular  momentum  due 
to  impulsive  forces  and  moments,  expressed  in  eqn.  (2.19).  The  impact  response  depended  upon 
impact  location  and  mass  of  the  impacting  hammer.  Generalized  displacement  results  are 
presented  for  several  impact  points  on  the  plate  using  and  impacting  hammer  of  about  10  lb. 
These  results  were  generated  for  a  plate  with  unequal  spring  stiffnesses,  discrete  masses  arranged 
arbitrarily,  and  damping,  which  is  case  IM-K5-M6-L4 

To  evaluate  the  initial  conditions  required  the  duration  of  the  impact  event  and  its 
intensity.  For  a  constant  impact  force,  eqn.  (2.19)  evaluates  to 

v((T)  =  ^ 
m 


0,(O+) 


MXXA  t 


(5.1) 


^(0+) 


M  ..At 


where  Fave  is  the  average  impact  force  intensity,  and  At  is  the  time  of  the  impact.  For  this 
analysis  the  impact  intensity  of  10  lb  for  a  duration  of  0.02  secs  was  used.  Impact  analysis  does 
not  necessary  start  at  time  t  =  0  secs.  With  the  use  of  the  kronecker  delta  function  5,  this  delay 
was  readily  incorporated.  Table  5.1  provided  locations  where  the  impact  would  take  place. 
These  locations  were  not  unique  and  represented  any  number  of  combinations  of  coordinates 
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where  an  impact  could  take  place.  Impact  locations  along  the  z-axis  produced  target  plots  that 
were  essentially  vertical.  The  nomenclature  introduced  in  chapter  three  for  tracking  vibration 
cases  is  augmented  here  to  track  impact  location.  The  new  notation  introduces  P;.  Pi  in  the  case 
IM-K5-M6-L4-P1  specifies  the  impact  point  to  correspond  to  that  presented  in  table  5.1.  All 
parameters  presented  for  damped  analysis  are  maintained  here. 


Table  5.1  Impact  Location  Placement 

Location 

(Xl,Zi) 

Pi 

(0,0) 

P2 

(18,0) 

P3 

(-18,0) 

The  case  considered  for  impact  analysis  was  IM-Ks-Me-LrPi .  The  initial  conditions  were 
vo(0+)  =  -0.0382848  ft/s,  0X(O+)  =  0.0  rad/s  and  0Z(O+)  =  0.0753376  rad/s.  The  solution  of 


eqn.(2.17)  was 
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u 0 (t)  —  (7.6899  xlO'7  e"1'13048  ("a5+t)  Sin[53.220  4  (-0.5  +  t)] 

+  2.4318  x  10"4e"a39845  (  a5+t)  Sin[31.723  5  (-0.5  +  t)] 

+  7.6868  x  10"5  e"0'399138  ("a5+t)  Sin[3 1.750  7  (-0.5  +  t)] 

+  1.6105  x  10"6  e1 09309  <  a5+t)Sin[52.341  9  (-0.5  +  t)] 

+  8.2506  xlO"8  e"1  12574  (a5+t)  Sin[53.109  8  (-0.5  +  t)])J(t -0.5) 

v 0 (t)  =  (-9.3065  xlO  8  e"1  13048  ("a5+t)  Sin[53.220  4  (-0.5  +  t)] 

-8.4358  x  10"5e"°'39845  ("a5+t)  Sin[31.723  5  (-0.5  +  t)] 

-1.5038  x  1 0  5  e-°-399138  (-0-5+0  Sin[3 1.750  7  (-0.5  +  t)] 

-5.5969  xlO  7  e  1  09309  (a5+t) Sin[52.341  9  (-0.5  +  t)] 

-2.8323  xlO  9  e  1  12574  <a5+t)  Sin[53.109  8  (-0.5  +  t)])£(t -0.5) 

w 0 (t)  —  (6.8405  xlO"7  e0, 13048  <a5+t)  Sin[53.220  4  (-0.5  +  t)] 

-1.8476  xlO  4 e  a39845  <a5+t)  Sin[31.723  5  (-0.5  +  t)] 

+  1.0800  xlO"4  e"°-399138  ("a5+t) Sin[3 1.750  7  (-0.5  +  t)] 

-1.5808  xlO"6  e"1  09309  (a5+t)Sin[52.341  9  (-0.5  +  t)] 

+  1.3300  xlO"7  e"1 12574  ("a5+t)  Sin[53.109  8  (-0.5  +  t)])£(t  -0.5) 


(5.2) 
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0x(t)  =  (-6.476 lx  10  6  e4'13048  <  a5+t)  Sin[53.220  4  (-0.5  +  t)] 

-1.9894  x  10  5  e  4)39845  (4)'5+t)  Sin[3 1.723  5  (-0.5  +  t)] 

+  9.9589  xlO'6  6  0399138  (-°5+t)  Sin[3 1.750  7  (-0.5  +  t)] 

+  1.3363  x  10"5  e4  09309  (4)'5+t)  Sin[52.341  9  (-0.5  +  t)] 

- 7.51  lx  iQ-7  e  4  12574  (“a5+t)  Sin[53.109  8  (-0.5  +  t)])£(t -0.5) 

6»y(t)  =  (-2.6473  x  10"6  e  4  13048  (4,'5+t)  Sin[53.220  4  (-0.5  +  t)] 

+  1.7463  x  1 0  "6  e "039845  (-°-5+t)  Sin[3 1.723  5  (-0.5  +  t)]  (5  3) 

-1.7869  x  10  6  64)399138  (  a5+t)Sin[3 1.750  7  (-0.5  +  t)] 

+ 1.005 lx  10‘6  e4  09309  (-°-5+t)  Sin[52.341  9  (-0.5  +  t)] 

+  1.6874  x  10  6  e4  12574  (a5+t)  Sin[53. 109  8  (-0.5  +  t)])£(t  -  0.5) 

0z(t)  =  (6.9201x  10  6  e  4  13048  (4)'5+t)  Sin[53.220  4  (-0.5  +  t)] 

-  2.6982  x  10"5e  4)39845  <  a5+t)  Sin[31.723  5  (-0.5  +  t)] 

-  7.3369  x  1 0 45  64,399138  (a5+t)Sin[3 1.750  7  (-0.5  +  t)] 

+  1.3258  x  10  5  e4  09309  (4X5+t)  Sin[52.341  9  (-0.5  +  t)] 

+  5.0184  x  10  7  e4  42574  (4)'5+t)  Sin[53.109  8  (-0.5  +  t)])£(t  -0.5) 

with  the  natural  frequencies  of  coni  =  5.0807  Hz,  con2=  5.0851  Hz,  con3=  5.1117  Hz,  con4=  8.4179 

Hz,  con5  =  8.5427  Hz,  con6=  8.5607  Hz,  and  the  damped  natural  frequencies  of  codi  =  5.0489  Hz, 

C0d2=  5.0532  Hz,  C0d3=  5.0847  Hz,  C0d4=  8.3304  Hz,  C0d5=  8.4526  Hz,  and  C0d6=  8.4702  Hz. 
Figure  5.1  presents  plots  of  displacement  vs.  time.  In  these  graphs,  displacements  Vo  and  Wo 
moved  in-phase  while  angular  rotations  0X  and  0Z  motions  were  in-phase  with  the  magnitude  of 
0y.  However,  0y  was  two  orders  of  magnitude  smaller  than  either  0X  or  0Z.  Figure  5.2  provides 
the  target  plot  for  this  case  showing  the  effect  of  angular  rotations  0X  and  0Z  on  this  traced 
displacement.  Figure  5.3  and  5.4  and  figures  5.5  and  5.6  are  displacement-target  plot  sets  for 


impact  locations  at  P2  and  P3,  respectively. 
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Figure  5.1  Displacement  plots  for  impact  unequal  stiffnesses,  added  discrete  with  a  =  0.001562, 

(3  =  0.001562 
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Figure  5.2  Target  plot  for  impact  analysis  with  impact  located  at  the  origin. 
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Figure  5.3  Displacement  plots  for  impact  at  location  IF 
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Figure  5.4  Target  plot  for  impact  at  location  Pr 
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Figure  5.5  Displacement  plots  for  impact  at  location  P3 
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Figure  5.6  Target  plot  for  impact  at  location  P3 
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CHAPTER  6 
FORCED  VIBRATION 

6.1  Introduction 

In  this  chapter,  forced  vibration  analysis  of  the  platform  system  is  presented.  Unlike 
previous  results,  the  forcing  function  was  constantly  applied  to  the  system,  which  produced  both 
a  transient  and  steady  state  solution.  Also  unlike  the  previous  discussions,  the  governing 
equations  of  motion  included  representations  of  the  forcing  functions,  depicted  as  inertial 
actuators.  This  inclusion  expanded  the  system  of  governing  equations  from  six  degrees  to  as 
many  as  ten  degrees  of  freedom.  The  equations  of  motion  were  solved  for  the  steady-state 
solutions  since  previous  results  captured  transient  results. 

Two  inertial  actuators  are  typical  for  the  optical  platform,  and  each  is  modeled  as  a 
second  order  system.  These  actuators  were  placed  and  oriented  in  such  a  way  to  generate  in¬ 
plane  displacement  uo,  transverse  displacement  vo,  and  all  three  angular  rotations.  The  forcing 
function  was  represented  as  a  single  harmonic  function  such  as  Fl  (t)  =  FlSin(colt)  or  a  linear 
combination  of  multiple  harmonic  functions  such  as 

F(t )  =  FxSin((0 \  t)  +  F2Sin{co2  t)  +  F2Sin(Q)2  t) .  Both  cases  were  considered  here.  Since  the 
problem  was  linear  in  nature,  displacement  solutions  for  the  case  of  multiple  frequency  forcing 
function  were  obtained  through  a  superposition  of  the  displacement  solution  at  each  forcing 
frequency.  Results  presented  here  considered  the  case  of  unequal  stiffenesses  with  added 
discrete  masses  arbitrarily  located  on  the  plate. 

Refer  to  chapter  4  on  free  vibration  for  an  explanation  of  the  systematic  notation  used 
within  this  chapter  to  identify  which  stiffness,  mass  composition,  and  mass  location  cases  are 
used  in  each  study.  These  variables  follow  the  pattern  outlined  in  tables  4.1  through  4.4. 
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6.2  Forced  Vibration  Analysis 

Forced  vibration  of  the  platform  used  two  actuators  whose  forcing  functions  were  given 
by 

J  Fj(0  =  10  Sin(J5t) 

F2(t)  =  10  Sin(25t)  (6.1) 

Forcing  function  Fi(t)  was  oriented  45°  relative  to  the  horizontal,  and  F2(t)  was  oriented  in  the 

vertical  position  only.  The  equations  that  governed  steady-state  forced  vibration  were  given  by 

(2.17).  These  equations  are  shown  in  complete  detail  in  Appendix  A-l.  A  modal  solution 

approach  was  used  to  decouple  the  problem,  solve  it  in  modal  coordinates,  and  transform  it  back 

to  the  original  coordinate  system  to  obtain  the  solution  of  the  original  problem.  The  system 

investigated  was  FF-K5-M6-L4.  The  constants  a  and  [3  used  for  the  damping  ratio  were 

0.0007812.  With  forcing  functions  given  by  eqn.  (6.1)  and  system  properties  given,  the  steady 

state  solutions  consisted  of  natural  frequencies  of  coni  =  4.9414  Hz,  con9  =  6.4366  Hz,  con3  = 

7.2966  Hz,  con4=  8.5514  Hz,  con5=  8.6176  Hz,  con6=  9.2568  Hz  and  damped  natural  of  C0di  = 

4.9114  Hz,  C0d2=  6.3855  Hz,  C0d3=  7.2310,  C0d4=  8.4612  Hz,  C0d5=  8.5259  Hz,  C0d6=  9.1511  Hz. 

The  displacement  equations  became 
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u0(t)  =  1.00406xl0"6  Sin[0.0121293  -75  t]  -  0.000321209  Sin[0.0240309  -75  t] 

+  1.31285xl0"6  Sin[0.0349573  -75  t]-5.42389xl0'7  Sin[0.0617196  -75  t] 

-9.55247xl0"6  Sin[0.0637154  -75  t]+5.38026xl0'6  Sin[0.0881979  -75  t] 

+4.3727xl0~6  Sin[0. 0198073  +25  t]-5.81755xl0'7  Sin[0.0198115  +25  t]  (6.2) 

-5.26029xl0"6  SinfO. 0198166  +25  t]-1.62503xl0'7  Sin[0.0239584  +25  t] 

-0.0000139863  Sin[0.0248248  +25  t]-2.58145xl0-6  Sin[0.0249298  +25  t] 

+6.37135xlO"7  Sin[0.0278009  +25  t]+0.0000994562  Sin[0.0316172  +25  t] 

-0.0000191568  SinfO.0555397  +25  t]+2.53003xl0-6  Sin[0.06701 12  +75  t] 

-0.0000395334  Sin[0.0671549  +75  t]+3.24891xl0-6  Sin[0.0673307  +75  t] 
v0(t)  =  0.000279697  Sin[0.0121293  -75  t]+1.00761xl0-6  Sin[0.0240309  -75  t] 

-1.94145xl0"6  SinfO. 0349573  -75  t]-2.69454xl0-7  Sin[0.0617196  -75  t] 

-8.67708xl0"7  SinfO. 0637154  -75  t]-8.28308xl0-6  SinfO.0881979  -75  t] 

+0.0000409697  SinfO.0198073  +25  t]+1.43775xl0-7  SinfO.01981 15  +25  t]  (6.3) 

+5.30289xl0"7  SinfO.0198166  +25  t]+2.50179xl0'7  SinfO.0239584  +25  t] 

-1.27046xl0"6  SinfO. 0248248  +25  t]-1.28244xl0-6  SinfO.0249298  +25  t] 

-9.42197xl0"7  SinfO. 0278009  +25  t]-3.1 1986xlO-7  SinfO.0316172  +25  t] 

-0.00533642  SinfO.0555397  +25  t]+0.0000237049  SinfO.0670112  +75  t] 

+9.77027xl0"6  SinfO.0671549  +75  t]-3.27522xl0-7  Sin[0.0673307  +75  t] 
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w0(t)=- 1 .444 1 7x  1 0"6  Sin[0.0121293  -75  t]-1.82636xl0-6  Sin[0.0240309  -75  t] 

-0.0000923278  Sin[0.0349573  -75  t]-2.91666xl0’7  Sin[0.0617196  -75  t] 

+0.0000194523  Sin[0.0637154  -75  t]+0.0001 15661  Sin[0.0881979  -75  t] 

-4.19896xl0"7  SinfO. 0198073  +25  t]-2.37918xl0'9  Sin[0.0198115  +25  t] 

+5.15085xl0"7  SinfO.0198166  +25  t]-3.49338xl0-6  Sin[0.0239584  +25  t]  (6.4) 

+0.000028481  SinfO.0248248  +25  t]-1.38816xl0-6  Sin[0.0249298  +25  t] 

-0.0000448072  Sin[0.0278009  +25  t]+5.65496xl0-7  Sin[0.0316172  +25t] 

+0.0000275536  Sin[0.0555397  +25  t]-2.4295xl0-7  Sin[0.06701 12  +75  t] 

-1.61678xl0"7  SinfO. 067 1549  +75  t]-3.18132xlO-7  Sin[0.0673307  +75  t] 

0x(t)=  4.60572xl0"6  SinfO.0121293  -75  t]-2.07747xl0-6  Sin[0.0240309  -75  t] 

+0.0000448641  SinfO.0349573  -75  t]+7.71 136xl0-6  SinfO.0617196  -75  t] 

+8.16241xl0"6  Sin[0.0637154  -75  t]+0.000256036  Sin[0.0881979  -75  t] 

+0.0000109739  SinfO.0198073  +25  t]+8.90035xl0-8  SinfO.01981 15  +25  t]  (6.5) 
-0.000013611  SinfO.0198166  +25  t]-7.7332xl0'6  SinfO.0239584  +25  t] 

+0.000011951  SinfO.0248248  +25  t]+0.0000367015  Sin[0.0249298  +25t] 

+0.0000217728  Sin[0.0278009  +25  t]+6.43246xl0-7  SinfO.0316172  +25  t] 

-0.0000878739  Sin[0.0555397  +25  t]+6.34946xl0-6  Sin[0.06701 12  +75  t] 

+6.04827xl0"6  Sin[0.0671549  +75  t]+8.40657xl0-6  SinfO.0673307  +75  t] 
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0y(t)=-2.82214xlO"8  Sin[0.0121293  -75  t]+6.7708xl0'6  Sin[0.0240309  -75  t] 

+0.0000196566  Sin[0.0349573  -75  t]-8.9916xl0'6  Sin[0.0617196  -75  t] 

+0.0000263297  Sin[0.0637154  -75  t]+0.0001 11047  Sin[0.0881979  -75  t] 

-1.41948xl0"6  SinfO. 0198073  +25  t]+1.76317xl0'7  Sin[0.0198115  +25  t] 

+  1.40466xl0"6  SinfO.0198166  +25  t]-3.35401xl0-6  Sin[0.0239584  +25  t]  (6.6) 

+0.0000385506  Sin[0.0248248  +25  t] -0.0000427947  Sin[0.0249298  +25  t] 

+9.53946xl0"6  Sin[0.0278009  +25  t]-2.09645xl0"6  Sin[0.0316172  +25  t] 

+5.38444xl0"7  Sin[0.0555397  +25  t]-8.21303xl0-7  Sin[0.06701 12  +75  t] 

+0.0000119817  SinfO.0671549  +75  t]-8.67559xlO-7  Sin[0.0673307  +75  t] 

0z(t)=  4.69373xl0"6  Sin[0.0121293  -75  t]-0.0000375903  Sin[0.0240309  -75  t] 

+4.25126xl0"6  Sin[0.0349573  -75  t]+1.42329xl0-6  Sin[0.0617196  -75  t] 

+0.000331645  Sin[0.0637154  -75  t] -0.0000274044  Sin[0.0881979  -75  t] 

-0.0000150431  Sin[0.0198073  +25  t]+3.08941xl0-7  Sin[0.01981 15  +25  t] 

-0.0000113176  SinfO.0198166  +25  t]+8.27713xl0-7  Sin[0.0239584  +25  t]  (6.7) 

+0.000485579  Sin[0.0248248  +25  t]+6.77404xl0-6  Sin[0.0249298  +25  tl 
+2.063 16xl0"6  Sin[0.0278009  +25  t]+0.00001 16391  Sin[0.0316172  +25  t] 

-0.000089553  Sin[0.0555397  +25  t]-8.70387xl0-6  Sin[0.06701 12  +75  t] 

+0.0000209942  Sin[0.0671549  +75  t]+6.9901xl0-6  Sin[0.0673307  +75  tl 
Although  long  and  tedious,  the  solutions  contained  no  more  than  seven  combinations  of  the  Sine 
function  at  each  frequency  and  a  phase  angle.  These  equations  were  plotted  in  figure  6.1.  Most 
notable  was  the  slowing  period  of  oscillation  in  displacement  vo(t)  and  rapid  oscillation  in  0z(t). 
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Figure  6.2  shows  the  same  input  information  with  the  exception  that  the  frequencies  were 
exchanged.  With  these  changes  it  was  possible  to  assess  the  effect  of  the  actuator  orientation  on 
the  motion  of  the  platform.  By  comparing  Figures  6.1  and  6.2  it  was  seen  that  the  actuator 
oriented  at  45°  dominates  the  motion  when  forcing  amplitudes  were  equal.  Figure  6.3  shows  the 
displacement  results  using  same  frequency  of  25  rad/s.  In-plane  displacement  Wo  was  relatively 
small  when  compared  to  the  other  two  components.  Also,  rotations  0X  and  0y  were  in-phase. 
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Figure  6.1  Displacement  plots  with  Fi(t)  =  10  Sin(75t)  at  45° 
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Figure  6.2  Displacement  plots  with  Fi(t)  =  10  Sin(25t)  at  45° 
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Figure  6.3  Displacement  plots  with  Fi(t)  =  F2(t)  =10  Sin(25t)  at  45° 


By  considering  that  the  forcing  functions  were  linear  combinations  of  harmonic  functions,  given 
by 

N 

Fk(t)  =  YJAklSin[(Oit ]  (6  8) 


it  was  possible  to  define 
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N 

Fl(t)  =  YJ  AuSin[Ci).t] 

i=  1 
N 

F2(t)  =  ^A,iSin\(d;t\ 

1=1 

Here  N  is  taken  as  7,  all  of  the  magnitudes  Aki  equal,  and  random  choices  of  frequency  were 
taken,  which  produced  the  following  forcing  functions 

Fx(t)  =  10[5m(5.99t)  +  5m(12.08t)  +  5m(16.45t)  +  Sin(35.0t)  +  Sin(43.\2t)  +  Sin(47.0St )  +  Sin(56.23t)\ 

Ft  (t)  =  10[5m(l  1.090  +  Sin(2\.llt)  +  5m(24.880  +  Sin(32.16t)  +  5m(41.990  +  Sin(41 .93t)  +  Sin(6().22t)\ 

Figures  6.4  and  6.5  show  time  history  plots  of  forcing  functions  Fi  and  F2. 
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(6.9) 


Figure  6.5  Plot  of  forcing  function  F2(t) 
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The  solution  process  for  this  case  merely  required  a  superimposed  solution  in  modal  space  and  a 


modal  transformation  back  to  the  original  coordinate  system.  The  results  were  quite  lengthy  and 
shown  only  graphically  below  in  figure  6.6 


VO 
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Figure  6.6  Displacement  plot  for  multiple  frequency  forcing  functions 
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CHAPTER  7 

CONCLUSION  AND  FUTURE  WORK 

This  project  studied  the  analysis  of  an  optical  platform  modeled  as  an  elastically  point 
supported  plate  containing  discrete  mass.  The  study  consisted  of  free  vibration,  impact  analysis 
and  force  vibration.  Both  undamped  and  damped  cases  were  considered,  where  the  damped 
cases  were  modeled  using  proportional  damping. 

1.  Dynamic  vibration  equations  of  equilibrium  were  derived  using  Lagrange’s  Equation. 

2.  Undamped  and  damped  analysis  was  considered 

3.  Two  types  of  motion  can  result  -  Direct  motion  due  to  initial  conditions  and  induced 
motion  due  to  system  characteristics.  Induced  motion  resembles  vibrational  beating 
-  slow  amplitude  changes  with  rapid  oscillations. 

4.  Equal  spring  stiffness  results  in  harmonic  motion,  without  beating  phenomena 

5.  Unequal  spring  stiffnesses  induce  motion  in  other  directions 

6.  Unequal  spring  stiffnesses  separate  modes  of  vibration 

7.  Added  masses  reduce  frequencies 

8.  Target  plot  density  is  greater  with  equal  spring  stiff  nesses  than  when  unequal. 

9.  Explicit  expressions  for  all  kinematic  variables  are  known,  leading  to  computation  of 
velocities  (linear  and  angular),  acceleartions  (linear  and  angular)  at  any  point  on  the 
plate. 

10.  All  experimental  results  can  be  compared  with  analytical  results 

There  are  several  ideas  that  can  be  considered  to  extend  the  work  of  this  project. 

1 .  Compare  analytical  results  with  experimental  results 
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2.  Understand  strategy  to  distribute  force  intensity  for  multiple  frequency  case 

3.  Develop  state  space  representation  for  future  control  applications 

4.  Incorporate  modeling  of  devices  such  as  fast  steering  mirror,  etc. 
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APPENDIX  A1 
EXPLICIT  EQUATIONS 


A  1.1  Coupled  Dynamic  Equation  for  Inertial  Actuators 
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A1.2  Coupled  Dynamic  Equations  for  Optical  Platform 
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APPENDIX  A2 
SUB  MATRICIES 

A2.1  Introduction 

This  appendix  presents  explicit  forms  of  the  matrices  used.  Recall  that  Lagrange’s 
Equation  provides  a  means  of  obtaining  a  system’s  governing  equations  of  motion,  which  in 
matrix  form  become 

[m]{x}  +  [c]{i} +  [&]{;*:}  =  {F}  (A2.1) 

Each  matrix  in  eqn.  A2.1  is  subdivided  into  four  matrices, 
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_[*>](6X4)  [^p  ]  (6x6)  _  |hj  1(6x1)  J  [  (6x1)  J 

where  [mj(4x4)  is  the  mass  matrix  of  the  shakers,  [mp](6x6)  is  the  mass  matrix  for  the  plate,  and 
its  discrete  masses,  [cj(4x4)  is  the  damping  matrix  of  the  shakers,  [cpJ(4x6)  is  the  damping  matrix 
that  couples  the  motion  of  the  plate  and  the  shakers,  [fcs](4v4)  is  the  stiffness  matrix  of  the  shakers, 


[&  ](4x6)is  the  stiffness  matrix  which  couples  the  motion  of  the  plate  and  the  shakers,  [k ](6x6)is 


the  stiffness  matrix  of  the  plate,  \  i<4  l)  l  is  a  vector  of  the  shaker  and  plate  accelerations, 

IMfixiJ 


\  4ll)l  is  a  vector  of  the  shaker  and  plate  velocities,  I  j  1  ,(4-'nl  js  a  vector  of  the  shaker  and 

P  1(6x1)  J  1  Wp  1(6x1)  J 


plate  displacements,  and  |  ^  j  is  a  vector  comprised  of  the  forces  imparted  by  the  shakers 


onto  the  system. 
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The  above  defined  matrices  for  the  shaker  are 
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Here  s  1  and  s2  stand  for  shaker  1  and  2  and  t  represents  the  portion  of  the  mass  of  the  shaker  in 
motion.  The  matrices  that  couple  the  motion  of  the  shaker  and  the  plate  are 
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The  matrices  that  govern  the  motion  of  the  plate  are 
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APPENDIX  B1 

MATHEMATICA  CODE  FOR  FORCED  VIBRATION  ANALYSIS 


ClearAll ; 

(* input*) 
f six  =  7.071; 
f sly=-7 . 071 ; 
f s2y=-10 ; 
tf=l .5; 

(Oil  =  5.99; 

C 012  =  12.08; 

0)13  =  16.45; 

(014  =  35.00; 

(015  =  43.12; 

(016  =  47.08; 

(017  =  75.0; 

(021=  11.09; 

(022  =  21.77; 

(023  =  24.88; 

(024  =  32.76; 

(025  =  41.99; 

(026  =  47.93; 

(027  =  60.22; 
(♦Stiffness  Terms*) 
a=0. 0007812; 

|3=0. 0007812; 
kl=115*12 ; 

(*kl=1612 .5 
k2=1612 . 5 
k3=1612 . 5 
k4=1612 . 5*) 
k3=120*12 ; 
k2=110*12 ; 
k4=105*12 ; 
(*kl=105*12 ; 
k2=125*12 ; 
k3=110*12 ; 
k4=115*12 ; *) 
ksl  =  3689.27; 
ks2=3689 . 27 ; 
mp=4 . 658 ; 
ma=0.155; (*ml*) 
mj=0.078; (*m2*) 
ml=  0.100;  ( *m3  * ) 
ml2=0 . 233; (*m4*) 
mtsl=  0.086335; 
mts2=  0.086335; 
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mbsl=  0.098447; 
mbs2=  0.098447; 

(*mp=l. 310488015;*) 

mt ot al=ma+m j  +ml+ml2  +mt  s 1 +mt  s  2  +mbs 1 +mbs  2  +mp ; 
(*xl  =  0.0  ; 
zl  =  0.0; 

yl=  0.0; 

xl2  =  0.0; 
zl2  =  0.0; 
yl2  =  0.0; 
xa  =0.0; 
za  =  0.0; 

ya  =  0.0; 

xj  =0.0; 
z j  =  0.0; 

yj  =  0.0;*) 

xl  =  -15/12;  (*  (6/7) * (1/2) *lplate; *) 

zl  =  -15/12;  (*  (6/7) *  (1/2) *wplate; *) 

yl=  (1/2) *hplate+0.333; 
xl2  =  6/12;  (* (6/7) * (1/2) *lplate; *) 

zl2  =  -6/12;  (*  (6/7) * (1/2) *wplate; *) 

yl2  =  (1/2) *hplate+0.333; 
xa  =9/12; 
za  =  9/12; 

ya  =  (1/2) *hplate+0 . 333; 
xj  =-12/12; 
zj  =  12/12; 

yj  =  (1/2) *hplate+0 . 333; 

xsl  =0;  (*  (6/7) * (1/2) *lplate; *) 

zsl  =- (2/7) * (1/2) *wplate; 
ysl  = (1/2) *hplate+0 . 333; 

xs2  =  -  (2/7) * (1/2) *lplate; 
zs2  =0 

ys2  = (1/2) *hplate+0 . 333; 

xbsl  =0;  (* (6/7) * (1/2) *lplate; *) 

zbsl  =- (2/7) * (1/2) *wplate; 
ybsl  = (1/2) *hplate; 

xbs2  =- (2/7) * (1/2) *lplate; 
zbs2  =0 

ybs2  = (1/2) *hplate; 

lplate=36/12; 
wplate=36/12; 
hplate  =  2/12; 

Ixx  = (mp/12) * (lplateA2+hplateA2)  ; 

Izz  =  (mp/12) * (wplateA2+hplateA2) ; 

Iyy  =  (mp/12) * (lplateA2+wplateA2) ; 

Ixz  =  0 ; 

Ixy  =  0; 
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lyz  =  0; 

zbar  =  0; 

xbar  =  0; 

ybar  =  0; 

lspring=17 . 25/12 ; 

wspring  =  17.25/12; 

hspring  =  (1/2) *hplate; 

zl  =  lspring; 

yl  =  -hspring; 

xl  =  wspring; 

z2  =  lspring; 

y2  =  -hspring; 

x2  =  -wspring; 

z3  =  -lspring; 

y3  =  -hspring; 

x3  =  -wspring; 

z4  =  -lspring; 

y4  =  -hspring; 

x4  =  wspring; 

position={usl [t] , vsl [t] , vs2 [t ] , uO [t] , vO [t] , wO [t] , 0x [t]  ,  0y  [t]  ,  0z  [t 

]}; 

velocity={usl ' [t] , vsl ' [t] , vs2 ' [t] , uO ' [t] , vO ' [t] , wO ' [t] , 0x ' [t] , 0y ' 
[t] ,0z' [t] }; 

acceleration^  usl ' '  [t]  ,  vsl ' '  [t]  ,  vs 2  '  '  [t  ]  ,  uO  '  '  [t]  ,  vO  ' '  [t]  ,  wO  ' '  [t]  , 
0x ' '  [t]  ,0y"  [t]  ,0z'  '  [t]  }  ; 

(*force={  f  six* Sin  [G)l*t  ]  ,  f  sly*Sin  [G)l*t]  ,  f  s2y*Sin  [0)2*t]  ,  f  s2z*Sin  [00 
2*t] ,0,0, 0,0, 0,0};*) 

F={  f  six*  (Sin  [G)ll*t]  +Sin  [0)12 *t  ]  +Sin  [00l3*t]  +Sin  [00l4*t]  +Sin  [00l5*t  ]  + 
Sin  [00l6*t]  +Sin  [0)17*t  ]  )  ,  f  sly*  (Sin  [G)ll*t  ]  +Sin  [00l2*t  ]  +Sin  [0)13*t  ]  +Si 


n  [0)14 *t]  +Sin  [COl 5*t]  +Sin  [C0l6*t]  +Sin  [C0l7*t]  )  ,  f  s2y*  (Sin  [0)21*t]  +Sin  [ 


C022*t  ]  +Sin  [0)23*t]  +Sin  [0)24*t]  +Sin  [0)25*t]  +Sin  [0)26*t]  +Sin  [0)27*t]  )  ,  0 


,  0,0,0,  0,0} 

B=  {{1,0, 0,0, 0,0, 0,0,0},  {0,1, 0,0, 0,0, 0,0,0},  {0,0, 1,0, 0,0, 0,0,0},  {0 
,0,0, 0,0, 0,0, 0,0},  {0,0, 0,0, 0,0, 0,0,0},  {0,0, 0,0, 0,0, 0,0,0},  {0,0,0, 
0,0, 0,0, 0,0},  {0,0, 0,0, 0,0, 0,0,0},  {0,0, 0,0, 0,0, 0,0,0},  {0,0, 0,0, 0,0 
,0,0,0}}; 

KEl= (1/2) *mtsl* ( (us 1 ' [t] ) A2+ (vsl ' [t] ) A2) ; 

KE2= (1/2) *mts2* ( (vs2 ' [t] ) A2) ; 


(*  Kinetic  Energy  of  the  Moving  Masses  of  Actuators*) 

(*  Potential  Energy  of  the  Springs  Attached  to  Actuators*) 


PE1= (ksl/2) * ( (usl [t] - (uO [t] +zsl*0y [t ] -ysl*0z [t] ) ) A2+ (vsl [t ] - 


(vO  [t]  —  z s  1  *0x  [t]  +xsl*0z  [t]  )  )  A2+  (wO  [t]  -xsl*0y  [t]  +ysl*0x  [t]  )  A2)  ; 
PE2= (ks2/2) * ( (uO [t] +zs2*0y [t] -ys2*0z [t] ) A2+ (vs2 [t] - (vO [t] - 


zs2*0x [t] +xs2*0z [t] ) ) A2+( (wO [t] -xs2*0y [t ] +ys2*0x [t ] ) ) A2) ; 

(*  Potential  Energy  of  the  Springs  Attached  to  Actuators*) 
(*  Potential  Energy  of  the  Springs  Attached  to  Actuators*) 


RD1= (csl/2) * ( (usl ' [t] - (uO ' [t ] +zsl*0y ' [t]-yl*0z' [t] ) ) A2+ (vsl ' [t ] - 
(vO ' [t] -zsl*0x ' [t] +xsl*0z ' [t] ) ) A2+ (wO ' [t] - 


xsl*0y ' [t] +ysl*0x ' [t] ) A2) ; 
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RD2= (cs2/2) * ( (uO ' [t ] +zs2*0y ' [t] -ys2*0z ' [t] ) A2+ (vs2 ' [t] - (vO ' [t ] - 
zs2*0x ' [t ] +xs2*0z ' [t] ) ) A2+( (wO' [t]-xs2*0y' [t]+ys2*0x' [t]))A2); 

(*  Potential  Energy  of  the  Springs  Attached  to  Actuators*) 

KE= (ml/2) * ( (uO ' [t ] +zl*0y ' [t]-yl*0z' [t])A2+(vO' [t ] - 
zl*0x ' [t ] +xl*0z ' [t ] ) A2+ (wO ' [t ] -xl*0y ' [t]+yl*0x' [t])A2)+ 

(ml2/2) * ( (uO ' [t] +zl2*0y ' [t] -yl2*0z ' [t] ) A2+ (vO ' [t] - 
zl2*0x ' [t ] +xl2*0z ' [t] ) A2+ (wO ' [t ] -xl2*0y ' [t]+yl2*0x' [t] ) A2)+ 

(mj/2) * ( (uO ' [t] +z j*0y ' [t]-yj*0z' [t])A2+(vO' [t ] - 
zj*0x' [t]+xj*0z' [t ] ) A2+ (wO ' [t] -x j*0y ' [t]+yj*0x' [t])A2)+ 

(ma/2) * ( (uO ' [t]+za*0y' [t]-ya*0z' [t])A2+(vO' [t ] - 
za*0x ' [t] +xa*0z ' [t])A2+(wO' [t]-xa*0y' [t]+ya*0x' [t])A2)+ 

(mbsl/2) * ( (uO ' [t ] +zbsl*0y ' [t ] -ybsl*0z ' [t] ) A2+ (vO ' [t] - 
zbsl*0x ' [t] +xbsl*0z ' [t] ) A2+ (wO ' [t] -xbsl*0y ' [t ] +ybsl*0x ' [t] ) A2) + 
(mbs2/2) * ( (uO ' [t ] +zbs2*0y ' [t ] -ybs2*0z ' [t] ) A2+ (vO ' [ t ] - 
zbs2*0x ' [t] +xbs2*0z ' [t] ) A2+ (wO ' [t] -xbs2*0y ' [t ] +ybsl*0x ' [t] ) A2) + 

(mp/2 ) * ( (uO '  [t])A2+(vO'  [t])A2+(wO'  [t])A2) +mp*xbar* (vO  '  [t ] *0z '  [ t ]  — 
wO '  [ t ] *0y '  [ t ] ) +mp*zbar* (uO  '  [t ] *0y '  [t ] - 

vO '  [t] *0x'  [t] )  +  (1/2) * ( Iyy*  (0y '  (t ] ) A2  +  Ixx*  (0x'  [t] ) A2  +  Izz*  (0z '  [t ] ) A 
2 ) -Ixz  *0x '  [t] *0z '  [t] ; 

PE= (kl/2) * ( (uO [t] +zl*0y [t] -yl*0z [t] ) A2+ (vO [t] - 
zl*0x [t] +xl*0z [t ] ) A2+ (wO [t] - 

xl*0y[t]+yl*0x[t] ) A2)+(k2/2) * ( (uO [t] +z2*0y [t] -y2*0z [t] ) A2+(v0[t]- 
z2*0x [t] +x2*0z [t ] ) A2+ (wO [t] - 

x2*0y[t]+y2*0x[t] ) A2)+(k3/2) * ( (uO [t ] +z3*0y [t] -y3*0z [t] ) A2+(v0[t]- 
z3*0x [t] +x3*0z [t ] ) A2+ (wO [t] - 

x3*0y[t]+y3*0x[t] ) A2)+(k4/2) * ( (uO [t ] +z4*0y [t] -y4*0z [t] ) A2+(v0[t]- 
z4*0x [t] +x4*0z [t] ) A2+(w0 [t] -x4*0y [t] +y4*0x [t] ) A2) ; 

(*  Potential  Energy  of  the  Springs  Attached  to  Plate*) 

(*  Rayleigh  Dissipation  Energy  of  the  Dampers  Attached  to  Plate*) 
RD= (cl/2) *  (  (uO '  [t]+zl*0'y[t] -yl*0z '  [t])A2+(vO'  [t ] - 
z 1 *0x ' [t]+xl*0z' [t ] ) A2+ (wO ' [t ] - 

xl*0y ' [t ] +yl*0x ' [t] ) A2) + (c2/2) * ( (uO ' [t]+z2*0y' [t ] - 
y2*0z ' [t ] ) A2+ (vO ' [t ] -z2*0x ' [t]+x2*0z' [t])A2+(wO' [t ] - 
x2*0y ' [t ] +y2*0x ' [t ] ) A2) + (c3/2) * ( (uO ' [t]+z3*0y' [t ] - 
y3*0z ' [t] ) A2+ (vO ' [t ] -z3*0x ' [t] +x3*0z ' [t] ) A2+ (wO ' [t ] - 
x3*0y ' [t] +y3*0x ' [t ] ) A2) + (c4/2) * ( (uO ' [t]+z4*0y' [t ] - 
y4*0z ' [t ] ) A2+ (vO ' [t ] -z4*0x ' [t]+x4*0z' [t])A2+(wO' [t ] — 
x4*0y '  [t ] +y4*0x '  [t] ) A2)  ; 

eqn [ i _ ] : =D [D [ (KE1+KE2+KE) , velocity [ [i] ] ] , t] - 

D [ (KE1+KE2+KE) , position [ [i] ] ] +D [ (RD1+RD2+RD) , velocity [ [i] ] ] +D [ (PE 
1+PE2+PE) , position [ [i] ] ] ; 

St if fness=Table [Coefficient [eqn [ j] , position [[k]]],{k,9},{j,9}]; 
Mass=Table [Coefficient [eqn [ j] , acceleration [ [k] ] ] , {k, 9} , { j,  9} ] ; 
Damping=Table [a*Mass [ [ j , k] ] +P*Stif fness [[j,k]],{k,9},{j,9}]; 
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L=Chop [CholeskyDecomposition [Mass] ] ; 

MDl=Inverse [Transpose [L] ] .Mass . Inverse [L] ; 

KDl=Inverse [Transpose [L] ] . Stiffness . Inverse [L] ; 

CDl=Inverse [Transpose [L] ] . Damping . Inverse [L] ; 

FDl=Chop [Inverse [Transpose [L] ] . F ] ; 
vecs=Chop [Eigenvectors [KDl] ] ; 

P=Transpose [vecs]  ; 

Ptranspose  =  vecs; 

S=Inverse[L] .P; 

KD2  =Chop [Ptranspose . KDl . P]  ; 

MatrixForm[%] ; 

MD2  =Chop [Ptranspose .MDl . P] ; 

MatrixForm[%] ; 

CD 2  =Chop [Ptranspose . CDl . P] ; 

MatrixForm[%] ; 

FD2=Chop [ Simplify [Ptranspose . FDl ] ] 

MatrixForm[%] ; 

particularll=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [coll*t ]  ]  ,  { j ,  9}  ]  ; 
particularl2=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [col2*t ]  ]  ,  {  j ,  9}  ]  ; 
particularl3=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [0)13*t ]  ]  ,  {  j ,  9}  ]  ; 
particular  14=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [G)14*t]  ]  ,  {  j ,  9}  ]  ; 
particularl5=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [G)15*t  ]  ]  ,  {  j ,  9}  ]  ; 
particularl6=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [col6*t ]  ]  ,  {  j ,  9}  ]  ; 
particularl7=Table  [Coefficient  [FD2  [  [  j]  ]  ,  Sin  [G)17*t  ]  ]  ,  {  j ,  9}  ]  ; 
MatrixForm[%] ; 

particular21=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [co21*t ]  ]  ,  {  j ,  9}  ]  ; 
particular22=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [co22*t ]  ]  ,  {  j ,  9}  ]  ; 
particular23=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [co23*t ]  ]  ,  {  j ,  9}  ]  ; 
particular24=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [co24*t ]  ]  ,  {  j ,  9}  ]  ; 
particular25=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [C025*t ]  ]  ,  {  j ,  9}  ]  ; 
particular26=Table  [Coefficient  [FD2  [  [  j ]  ]  ,  Sin  [co26*t ]  ]  ,  {  j ,  9}  ]  ; 
particular27=Table  [Coefficient  [FD2  [  [  j]  ]  ,  Sin  [co27*t ]  ]  ,  {  j ,  9}  ]  ; 
MatrixForm[%] ; 

vals=Chop [Eigenvalues [KDl] ]  ; 

C0n=Sqrt  [vals]  ; 

C=(a/  (2*C0n)  )  +  (P*C0n)  /2; 

G)d=Table[G)n[  [  j]  ]  *Sqrt  [l-^[  [  j]  ]  ]  ,  { j,  9}]  ; 
fn=Sqrt [vals] / (2*Pi) ; 
fd=C0d/  (2*Pi)  ; 

{23.71629340343848'  (Sin[11.09'  t] +Sin [21 . 77'  t] +Sin [24 . 88' 

t] +Sin [32.76'  t] +Sin [41 . 99'  t ] +Sin [47 . 93'  t] +Sin [60 . 22' 

t] )+12. 176380963697648'  (Sin[5.99'  t] +Sin [12 . 08'  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t] +Sin [43 . 12 '  t ] +Sin [47 . 08'  t]+Sin[75.' 

t] ), 0.47429779935456384'  (Sin[11.09'  t] +Sin [21 . 77'  t] +Sin [24 . 88' 

t] +Sin [32 .76'  t] +Sin [41 . 99'  t ] +Sin [47 . 93'  t] +Sin [ 60 . 22' 

t] ) +28. 545273262255563'  (Sin[5.99'  t] +Sin [12 . 08'  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t] +Sin [43 . 12 '  t ] +Sin  [47 . 08'  t]+Sin[75.'  t]),- 

23.912619107296337'  (Sin[11.09'  t] +Sin [21 . 77'  t] +Sin [24 . 88' 

t] +Sin [32 .76'  t] +Sin [41 . 99'  t ] +Sin [47 . 93'  t] +Sin [60 . 22' 


t] )+13. 049440953020028'  (Sin[5.99'  t ] +Sin [12 . 08'  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t ] +Sin [47 . 08'  t]+Sin[75.'  t]),- 
0.055967219247049844'  (Sin[11.09'  t] +Sin [21 . 77'  t ] +Sin [24 . 88' 

t] +Sin [32 .76'  t ] +Sin [41 . 99'  t ] +Sin [47 . 93'  t ] +Sin [ 60 . 22 '  t])- 
1.511881440302727'  (Sin[5.99'  t] +Sin [12 . 08'  t] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t ] +Sin [47 . 08'  t]+Sin[75.' 
t] ), 2.2775481852129817'  (Sin[11.09'  t] +Sin [21 . 77'  t] +Sin [24 . 88' 

t] +Sin [32 .76'  t] +Sin [41 . 99'  t ] +Sin [47 . 93'  t ] +Sin [ 60 . 22 '  t])- 
1.8192557459641425'  (Sin[5.99'  t] +Sin [12 . 08'  t] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t] +Sin [47 . 08'  t]+Sin[75.' 
t] ), 0.31869739985449297'  (Sin[11.09'  t] +Sin [21 . 77'  t]+Sin[24.88 

t] +Sin [32 .76'  t ] +Sin [41 . 99'  t ] +Sin [47 . 93'  t] +Sin [ 60 . 22'  t])- 
0.08118442721821306'  (Sin[5.99'  t ] +Sin [ 12 . 08 '  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t] +Sin [47 . 08'  t]+Sin[75.'  t]),- 
0.1684652471422208'  (Sin[11.09'  t] +Sin [21 . 77'  t ] +Sin [24 . 88' 

t] +Sin [32 .76'  t ] +Sin [41 . 99'  t ] +Sin [47 . 93'  t ] +Sin [ 60 . 22 ' 
t] )+0. 828263591485338'  (Sin[5.99'  t] +Sin [12 . 08'  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t ] +Sin [47 . 08'  t]+Sin[75.' 
t] ), 0.23605585449666983'  (Sin[11.09'  t] +Sin [21 . 77'  t]+Sin[24.88 

t] +Sin [32 .76'  t ] +Sin [41 . 99'  t ] +Sin [47 . 93'  t ] +Sin [ 60 . 22 ' 
t] )+3. 008949946399901'  (Sin[5.99'  t] +Sin [12 . 08'  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t] +Sin [47 . 08'  t]+Sin[75.' 
t] ), 4 .289741825716503'  (Sin[11.09'  t] +Sin [21 . 77'  t] +Sin [24 . 88' 

t] +Sin [32 . 76'  t ] +Sin [41 . 99'  t ] +Sin [47 . 93'  t ] +Sin [ 60 . 22 ' 
t] )+3. 087057744212988'  (Sin[5.99'  t] +Sin [12 . 08'  t ] +Sin [16 . 45' 

t] +Sin [35 . '  t ] +Sin [43.12'  t ] +Sin [47 . 08'  t]+Sin[75.'  t] ) } 
(♦Displacement  Solution  -  Each  solution  will  be  a  superposition 
of  two  solutions  because  od  the  two  forcing  functions*) 
Xll=Table  [particularll  [  [n]  ]  /Sqrt  [  (con  [  [n]  ]  A2- 
©11A2 )  A2+  (2*£[  [n]  ]  *©n  [  [n]  ]  *©11)  A2]  ,  {n,  1,9}]; 

011=Table  [ArcTan  [  (2*£  [  [n]  ]  *©n  [  [n]  ]  *©11)  /  (con  [  [n]  ]  A2- 
©11A2) ] , {n, 1, 9} ] ; 

Xl2=Table  [particularl2  [  [n]  ]  / Sqrt  [  (COn  [  [n]  ]  A2- 
©12A2)  A2+  (2*£[  [n]  ]  *C0n  [  [n]  ]  *©12)  A2  ]  ,  {n,  1 ,  9 }  ]  ; 

012=Table  [ArcTan  [  (2*£ [  [n]  ]  *© n  [  [n]  ]  *©12)  /  (con  [  [n]  ]  A2- 
©12A2) ] , {n, 1, 9} ] ; 

Xl3=Table  [particularl3  [  [n]  ]  / Sqrt  [  (G)n  [  [n]  ]  A2- 
©13A2) A2+  (2*£[  [n]  ]  *©n[  [n]  ]  *©13)  A2]  ,  {n,  1,  9}  ]  ; 

013=Table  [ArcTan  [  (2*£ [  [n]  ]  *©n  [  [n]  ]  *©13)  /  (©n  [  [n]  ]  A2- 
©13A2) ] , {n, 1, 9} ] ; 

Xl4=Table [particularl4 [ [n] ] / Sqrt [ (©n [ [n] ] A2- 
©14A2) A2+(2*^[ [n]  ]*©n[ [n] ]*©14) A2]  ,  {n,l,  9}]  ; 

014=Table  [ArcTan  [  (2*£  [  [n]  ]  *©n  [  [n]  ]  *©14)  /  (©n  [  [n]  ]  A2- 
©14A2) ] , {n, 1, 9} ] ; 

Xl5=Table [particularl5 [ [n] ] /Sqrt [ (©n [ [n] ] A2- 
©15A2) A2+ (2*^[ [n] ] *©n [ [n] ] *©15) A2] , {n, 1, 9} ] ; 

015=Table  [ArcTan  [  (2*^  [  [n]  ]  *©n  [  [n]  ]  *©15)  /  (©n  [  [n]  ]  A2- 
©15A2) ] , {n, 1 , 9} ] ; 
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Xl6=Table  [particularl6  [  [n]  ]  / Sqrt  [  (con  [  [n]  ]  A2- 
C0l6A2)  A2+  (2*£[  [n]  ]  *©n[  [n]  ]  *©16)  A2 ]  ,  {n,  1 ,  9 }  ]  ; 

016=Table  [ArcTan  [  (2*£  [  [n]  ]  *G)n  [  [n]  ]  *©16)  /  (con  [  [n]  ]  A2- 
©16A2) ] , {n, 1, 9} ] ; 

Xl7=Table  [particularl7  [  [n]  ]  /Sqrt  [  (G)n  [  [n]  ]  A2- 
0>17A2)  a2+  (2*£[  [n]  ]  *<On[  [n]  ]  *©17)  A2]  ,  {n,  1,9}]; 

017=Table  [ArcTan  [  (2*£  [  [n]  ]  *©n  [  [n]  ]  *©17)  /  (©n  [  [n]  ]  A2- 
©17A2) ] , {n, 1, 9} ] ; 

X21=Table [particular21 [ [n] ] /Sqrt [ (©n [ [n] ] A2- 
©21A2) A2+ (2*£[ [n]  ]  *©n [ [n] ] *©21) A2]  ,  {n,  1,  9}  ]  ; 

021=Table  [ArcTan  [  (2*£  [  [n]  ]  *©n  [  [n]  ]  *©21)  /  (©n  [  [n]  ]  A2- 
©21A2 )],{n,l,9}]; 

X22=Table [particular22 [ [n] ] / Sqrt [ (©n [ [n] ] A2- 
©22A2 ) A2+  (2*£[ [n] ] *©n[ [n] ] *©22) A2] ,  {n, 1, 9} ] ; 

022=Table  [ArcTan  [  (2*£ [  [n]  ]  *©n  [  [n]  ]  *©22)  /  (©n  [  [n]  ]  A2- 
©22A2 )],{n,l,9}]; 

X23=Table [particular23 [ [n] ] / Sqrt [ (©n [ [n] ] A2- 
©23A2 ) A2+ (2*£[ [n]  ]  *©n[ [n] ] *©23) A2] , {n,  1,  9} ]  ; 

023=Table  [ArcTan  [  (2*£ [  [n]  ]  *©n  [  [n]  ]  *©23)  /  (©n  [  [n]  ]  A2- 
©23A2) ] , {n, 1, 9} ] ; 

X24=Table [particular24 [ [n] ] /Sqrt [ (©n [ [n] ] A2- 
©24A2) A2+ (2*£[ [n] ] *©n[ [n] ] *©24) A2] , {n, 1, 9} ] ; 

024=Table  [ArcTan  [  (2*£  [  [n]  ]  *©n  [  [n]  ]  *©24)  /  (©n  [  [n]  ]  A2- 
©24A2) ] , {n, 1, 9} ] ; 

X25=Table [particular25 [ [n] ] /Sqrt [ (©n [ [n] ] A2- 
©25A2) A2+(2*^[ [n] ]*©n[ [n] ]*©25) A2] , {n,l,9}] ; 

025=Table  [ArcTan  [  (2*£  [  [n]  ]  *©n  [  [n]  ]  *©25)  /  (©n  [  [n]  ]  A2- 
©25A2) ] , {n, 1, 9} ] ; 

X26=Table [particular26 [ [n] ] / Sqrt [ (©n [ [n] ] A2- 
©26A2) A2+ (2*£[ [n] ] *©n[ [n] ] *©26) A2] , {n, 1, 9} ] ; 

026=Table  [ArcTan  [  (2*£ [  [n]  ]  *©n  [  [n]  ]  *©26)  /  (©n  [  [n]  ]  A2- 
©26A2) ] , {n, 1, 9} ] ; 

X27=Table [particular27 [ [n] ] / Sqrt [ (©n [ [n] ] A2- 
©27A2) A2+ (2*^[ [n] ] *©n[ [n] ] *©27) A2] , {n, 1, 9} ] ; 

027=Table  [ArcTan  [  (2*^  [  [n]  ]  *©n  [  [n]  ]  *©27)  /  (©n  [  [n]  ]  A2- 
©27A2) ] , {n, 1 , 9 } ] ; 

rxp=Chop [Table [Xll [ [k] ] *Sin [©ll*t+011 [ [k] ] ] +X12 [ [k] ] *Sin [©12*t+0 
12 [  [k]  ]  ]  + 

X13  [  [k]  ]  *Sin  [©13*t+013  [  [k]  ]  ]  +X14  [  [k]  ]  *Sin  [©14*t+014  [  [k]  ]  ]  +X15  [  [k 
] ] *Sin [©15*t+015 [ [k] ] ] +X16 [ [k] ] *Sin [©16*t+016 [ [k] ] ] +X17 [ [k] ] *Sin 
[©17*t+017 [ [k] ] ]+ 

+X21 [ [k] ] *Sin [©21*t+021 [ [k] ] ] +X22 [ [k] ] *Sin [©22*t+022 [ [k] ] ] + 
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X23  [  [k]  ]  *Sin  [G)23*t+023  [  [k]  ]  ]  +X24  [  [k]  ]  *Sin  [G)24*t+024  [  [k]  ]  ]  +X25  [  [k 
]  ]  *Sin  [oo25*t+025  [  [k]  ]  ]  +X26  [  [k]  ]  *Sin  [O)26*t+026  [  [k]  ]  ]  +X27  [  [k]  ]  *Sin 
[©27*t+027 [ [k] ]  ]  ,  {k, 1, 9} ] ] ; 

Sinverse=Ptranspose . L; 


xx=Chop [S . rxp] ; 
usll=xx [ [1] ] ; 
vsll=xx[ [2] ] ; 
vs22=xx [ [3] ] ; 
uO=xx [ [4] ] ; 
vO=xx [ [5] ] ; 
wO=xx [ [6] ] ; 

(|)x=xx  [  [  7  ]  ]  ; 

(|)y=xx  [  [8]  ]  ; 

(|)z=xx[  [9]  ]  ; 

plotf slx=Plot  [F  [  [1]  ]  ,  {t,  0,  tf } ,  PlotStyle— *Thin,  PlotRange— *A11,  Plot 
Label— *Subscript  ["F"  ,  1]  ,  AxesLabel-*  { secs,  lb}  ]  ; 

plotf sly=Plot  [F  [  [3]  ]  ,  {t,  0,  tf } ,  PlotStyle— *Thin,  PlotRange— *A11,  Plot 

Label— ^Subscript  [  "F"  ,  2]  ,  AxesLabel-*  { secs,  lb}  ]  ; 

plotf s2y=Plot  [F  [  [3]  ]  ,  {t,  0,  tf } ,  PlotStyle— *Thin,  PlotRange— *A11] 

Show  [plotf  six,  plotf  sly,  plotf  s2y] 

plotuO=Plot  [uO ,  {t ,  0 ,  tf } ,  PlotStyle— *Thin,  PlotRange— *A11,  PlotLabel— * 

Subscript  [u,  0]  ,  AxesLabel— »{ secs,  ft}  ] 

plotvO=Plot  [vO ,  {t ,  0 ,  tf } ,  PlotStyle— *Thin,  PlotRange- 

>A11 ,  PlotLabel— ^Subscript  [v,  0]  ,  Axes  Label— »{ secs,  ft}  ] 

plotw0=Plot  [wO ,  {t ,  0 ,  tf } ,  PlotStyle— >Thin,  PlotRange—»All,  Plot  Label— > 

Subscript  [w,  0]  ,  AxesLabel— >{ secs,  ft}  ] 

plot(|)x=Plot  [(|)x,  {t,  0,  tf } ,  PlotStyle— >Thin,  PlotRange— >A11,  PlotLabel— > 
Subscript  [0,  x]  ,  AxesLabel-* { secs,  rads }  ]  ; 

plot(j)y=Plot  [(|)y,  {t,  0,  tf } ,  PlotStyle— »Thin,  PlotRange— >A11,  PlotLabel— > 
Subscript  [0,  y]  ,  AxesLabel— »{ secs,  rads }  ]  ; 

plot(j)z=Plot  [ <|)z ,  {t,  0,  tf } ,  PlotStyle— »Thin,  PlotRange— >A11,  PlotLabel—* 
Subscript  [0,  z]  ,  AxesLabel— »{ secs,  rads }  ]  ; 

Show  [Graphics  Array  [  {  { plot  uO ,  plot  (j)x  } ,  {plotvO,  plot(j)y }  ,  {plotwO,  plot(|) 
z  }  }  ]  ] 


